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ABSTRACT 

We present maps, source lists, and derived number counts from the largest, unbiassed, 
}~^ • extragalactic sub-mm survey so far undertaken with the SCUBA camera on the JCMT. 

Our maps are located in two regions of sky (ELAIS N2 and Lockman-Hole E) and cover 
260 arcmin^, to a typical rms noise level of crgso — 2.5 mJy/beam. We have reduced 
the data using both the standard JCMT SURF procedures, and our own IDL-based 
pipeline which produces zero-footprint maps and noise images. The uncorrelated noise 
maps produced by the latter approach have enabled us to apply a maximum-likelihood 

O [ method to measure the statistical significance of each peak in our maps, leading to 

properly-quantified errors on the flux density of all potential sources. 

^ ■ We detect 19 sources with S/N > 4, and 38 with S/N > 3.5. To assess both 

' the completeness of this survey, and the impact of source confusion as a function of 

flux density we have applied our source-extraction algorithm to a series of simulated 
images. The result is a new estimate of the sub-mm source counts over the flux-density 

^ . , range 5*850 — 5 — 15 mJy, which we compare with estimates derived by other workers, 

jrt ' and with the predictions of a number of models. 

Our best estimate of the cumulative source count at 5*850 > 8 mJy is 320tioo per 
square degree. Assuming that the majority of sources lie at z > 1.5, this result implies 
that the co-moving number density of high-redshift galaxies forming stars at a rate 
in excess of 1000 MQyr~^ is ~ 10~^ Mpc~^, with only a weak dependence on the 
precise redshift distribution. This number density corresponds to the number density 
of massive ellipticals with L > 3 — 4L* in the present-day universe, and is also the same 
as the co-moving number density of comparably massive, passively-evolving objects 
in the redshift band 1 < z < 2 inferred from recent surveys of extremely red objects. 
Thus the bright sub-mm sources uncovered by this survey can plausibly account for 
the formation of all present-day massive ellipticals. Improved redshift constraints, and 
ultimately an improved measure of sub-mm source clustering can refine or refute this 
picture. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: starburst - infrared: galaxies 
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1 INTRODUCTION 

Over the past three years deep sub-milhmetre surveys, car- 
ried out to a variety of depths (Smail et al. 1997, Hughes et 
al. 1998, Barger, Cowie & Sanders 1999, Blain et al. 1999, 
Eales et al. 2000, Borys et al. 2002) have highhghted the ex- 
treme importance of dust in the determination of the global 
star-formation history of the Universe. Optical/UV studies 
originally suggested a steep rise in the star-formation rate 
(SFR) as a function of redshift between 2 = and 2 = 1 
(Lilly et al. 1996), peaking at 1.0 < z < 1.5, and declin- 
ing to values comparable to those observed in the present 
day at 2 ~ 4 (Madau et al. 1996). More recent optical/UV 
analyses have implied a much gentler trend in SFR density 
with epoch, Cowie et al. (1999) finding a more gradual slope 
(~ {l-\-zY'^ for an Einstein de Sitter cosmology) at low red- 
shift {z < 1) than had previously been determined by Lilly 
et al. (1996), and at high redshift the Lyman break pop- 
ulation (Steidel et al. 1999) suggesting an approximately 
uniform value for 1 < 2 < 4. However, in heavily dust- 
enshrouded star-forming regions, much of the optical/UV 
radiation emitted by the young stars is absorbed, leading to 
the possibility that a significant amount of star formation, 
particularly at high redshift, may have been missed in these 
wavebands. 

In order to quantify the star-formation density con- 
tributed by such highly-obscured objects, we must observe 
directly the rest-frame far-infrared (FIR) emission from the 
heated dust, which at redshifts greater than z ~ 1 is shifted 
into the sub-millimetre waveband. The steep spectral index 
of the thermal emission longwards of the peak at A ~ 100 /im 
results in a sufficiently large negative K-correction that the 
dimming effect of increasing cosmological distance is effec- 
tively cancelled over a wide range in redshift. Consequently, 
for a galaxy with fixed intrinsic FIR luminosity, we would 
expect to observe approximately the same 850 /im flux den- 
sity for an object at z = 8 as at z = 1. The strength of this 
effect is obviously dependent on the assumed cosmology. It is 
most pronounced in an Einstein-de Sitter universe, in which 
the predicted flux density of an object of fixed luminosity 
actually increases slightly beyond z ~ 1. If instead we adopt 
fijM = 0.3, Oa = 0.7, a very gentle decline in flux density is 
predicted over this same redshift range. 

Indications that a large fraction of star-forming activ- 
ity at high redshift is heavily enshrouded in dust have been 
provided, for example, by a deep sub-millimetre survey, of 
the northern HDF (Hughes et al. 1998). The discovery of 

5 discrete sources with either very faint or presently unde- 
tected optical identifications, accounting for approximately 
30 — 50% of the sub-millimetre background observed by 
COBE-FIRAS (Puget et al. 1996, Fixsen et al. 1998, Hauser 
et al. 1998), have suggested a star- formation density at least 
a factor of ~ 5 higher at z > 2 than that originally de- 
duced from optical/UV data. The results of Barger, Cowie 

6 Richards (2000) and Eales et al. (2000) are also consistent 
with a high redshift SFR density up to an order of mag- 
nitude higher than that inferred in the optical/UV. These 
sources, are believed to be analogous to the local ultralumi- 
nous infrared galaxies (ULIRGs) and undergoing a period 
of massive star-formation at rates of 10^ — lO^Meyr^^. 

The discovery of such objects at high redshift has im- 
portant implications for the formation mechanism of mas- 



sive elliptical galaxies. In current cold dark matter (CDM) 
based models, elliptical galaxies are created at modest red- 
shift by the merging of intermediate-mass discs (Baugh, Cole 
& Frenk 1996). Although some recent semi-analytic results 
appear supportive of the formation of massive ellipticals at 
modest redshift by means of hierarchical mergers (Kauff- 
man & Chariot 1998), observations continue to imply that 
a massive coeval starburst of timescale ~ 1 Gyr at high 
redshift (2 > 3) is required to explain the properties of at 
least some ellipticals (Jimenez et al. 1999, Renzini & Cimatti 
2000). Additionally, the connection between black-hole and 
spheroid mass (Magorrian et al. 1998, Kormendy & Geb- 
hardt 2001) suggests that galaxies hosting an active nucleus 
may be more representative of the general population than 
had previously been thought, and might indicate a link be- 
tween spheroid formation and the epoch of maximum quasar 
activity at 2 = 2 — 3. 

The discovery of the Lyman-break population ai z — 
2 — 4 (Steidel et al. 1999) has prompted the suggestion that 
the formation of present-day galactic bulges and ellipticals 
has already been observed at optical wavelengths. However, 
the most massive ellipticals contain 10^^ solar masses of 
stars, requiring a sustained SFR of ~ 1000 M0yr~^ for for- 
mation over a 1 Gyr period. Even if the most optimistic 
corrections are made for the effects of dust, the inferred star 
formation rates in these objects (typically 3 — 30 MQyr~^) 
fall short by 1 — 2 orders of magnitude. Moreover, the 
first direct tests of such corrections suggest they may in- 
deed have been exaggerated. Direct measurements of the 
sub-millimetre/millimetre emission of the strongly-lensed 
Lyman-break galaxies MS 1512+36-cB58 (van der Werf et 
al. 2001, Baker et al. 2001, Sawicki 2000) and MS 1358-f62- 
Gl (van der Werf et al. 2001) imply that the procedure used 
to predict the sub-millimetre emission of the Lyman-break 
population over-predicts the observed 850 pm flux densi- 
ties by a factor of up to 14. A recent X-ray detection of 
MS 1512-I-36-CB58 (Almaini, Pettini & Steidel in prep) also 
supports the view that corrections for dust have been exag- 
gereated, the level of X-ray emission implying a much lower 
SFR than that derived from the "dust-corrected" UV spec- 
trum, but consistent with the measured 850 pim flux density, 
for a typical starburst SED. 

One possibility is that spheroid formation is too widely 
distributed in space and/or time to give rise to such spectac- 
ular starbursts. However, observationally we may be miss- 
ing massive starbursts in the optical as a result of dust ob- 
scuration or because they are at too great a redshift. Sub- 
millimetre detections of the 2 ~ 4 radio galaxies 4C41.17 
(Dunlop et al. 1994, Hughes, Dunlop & Rawlings 1997) and 
8C1435-I-635 (Ivison 1995, Ivison et al. 1998) confirm that 
the Lyman-break population does not tell the full story. 
These extreme radio galaxies display the properties expected 
of youthful massive ellipticals, and have inferred dust masses 
> 1O®M0 and SFRs > 1000 MQyr-\ 

Here we present the first results from the 'SCUBA 
8 mjy Survey', which is a wider-area, somewhat shallower 
survey than its earlier counterparts (HDF-N - Hughes et al. 
1998; Clusters Survey - Smail et al. 1997, Blain et al. 1998; 
Hawaii Fields - Barger, Cowie & Sanders 1999; CUDSS - 
Eales et al. 2000), undertaken with the aim of constraining 
the brighter end of the 850 /im source counts in the region 
of 8 mJy. This survey has three key advantages over pre- 
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vious sub-millimetre surveys. Firstly, the relatively bright 
flux-density limit facilitates follow-up with existing instru- 
mentation which should eventually yield an accurate astro- 
metric position for every source via either the VLA, or in the 
case of the very high redshift sources (z > 3 — 4) IRAM PdB. 
Secondly, the steep slope of the source counts makes confu- 
sion a serious issue when probing deeper than ~ 2 — 3 mjy. 
The number density of sources brighter than 8 mJy is a fac- 
tor of ~ 3 lower than at 3 mJy and so source-confusion is 
much less of a problem in this survey. Finally, any sources 
discovered in this survey must be as bright or brighter in the 
sub-millimetre than the extreme radio galaxies mentioned 
above. This survey is therefore optimally placed for detect- 
ing the most luminous starbursts with SFR ~ 1000 M0yr~^, 
and through an extensive follow up program already under- 
way (Fox et al. 2002, Lutz et al. 2001) will provide clues to 
the formation of the most massive ellipticals. 

This paper is laid out as follows. In Sections 2 and 3 
we describe the observations and data reduction. In Section 
4 we outline the source extraction algorithm used to create 
our source catalogues. This is described in full mathematical 
detail in Appendix A. In Section 5 we present and discuss 
the results of simulations carried out both in conjunction 
with the real data, and through the production of fully- 
simulated survey areas. In Section 6 we present our source 
lists, giving positions, along with 850 and 450 /im flux densi- 
ties. In Section 7 we derive the resulting source counts (both 
raw and corrected using the simulation results). In section 
8 we present 2-point auto-correlation functions for each of 
our survey fields. A full discussion of the results, including 
comparison with other surveys, calculated dust masses, star- 
formation rates, estimates of co-moving number density, and 
clustering properties is given in Section 9. Finally, our con- 
clusions are summarised in Section 10. In subsequent papers, 
the cross-correlation of our SCUBA sources with deep X-ray 
and radio observations will be discussed by Almaini et al. 
(2002) and Ivison et al. (2002) respectively. 

For the calculation of physical quantities we have as- 
sumed Hq = 67 kms^^Mpc^^ throughout, and all cosmolog- 
ical calculations have been performed for both an Einstein- 
de Sitter universe and a A-dominated flat cosmology with 
Om = 0.3, f^A = 0.7. 



2 OBSERVATIONS 

The survey is divided between two flelds; the Lockman 
Hole East (centred at RA 10:52:8.82, DEC 4-57:21:33.8) and 
ELAIS N2 (centred at RA 16:36:48.85, DEC +41:01:48.5). 
These areas both have a vast quantity of multi-wavelength 
data available for follow-up studies, and were selected to co- 
incide with deep Infrared Space Observatory (ISO) surveys 
at 6.7, 15, 90 and 175 /im (Lockman Hole East - Elbaz et 
al. 1999, Kawara et al. 1998; ELAIS N2 - Oliver et al. 2000, 
Serjeant et al. 2000, Efstathiou et al. 2000). 

Pilot observations of the Lockman Hole East and ELAIS 
N2 fields began in March 1998 and July 1998 respectively, 
using the Sub-millimetre Common User Bolometer Array 
(SCUBA) (Holland et al. 1999) on the James Clerk Maxwell 
Telescope (JCMT) on Mauna Kea, Hawaii. In brief, SCUBA 
consists of two feedhorn-coupled bolometer arrays, both of 
which are arranged in a close-packed hexagon and have ap- 



proximately a 2.3 arcmin field of view (slightly smaller on 
the short-wavelength array). The diffraction- limited beams 
delivered by the JCMT have FWHM of 14.5" at 850 ^m 
and 7.5" at 450 pm, and the bolometer feedhorns on the ar- 
rays are sized for optimal coupling to the respective beams. 
As a result, the long-wavelength (LW) array consists of 37 
bolometers, and the short-wavelength (SW) array consists 
of 91 bolometers. Observations may be made simultane- 
ously with both arrays by means of a dichroic beamsplit- 
ter. A ^He/'^He dilution refrigerator cools SCUBA down 
to ~ 90 mK, making observations with the array sky- 
background limited. 

The most efficient mapping technique for sources 
smaller than the field of view is 'jiggle-mapping', and is the 
method employed in this survey. As explained above, the 
bolometers are arranged such that the sky is under-sampled 
at any instant and so the secondary is 'jiggled' in a hexag- 
onal pattern to fill in the gaps. In order to fully sample the 
sky at both wavelengths a 64-point jiggle-map is required. 
Temporal variations and linear spatial gradients in the sky 
were removed by chopping and nodding. A 30-arcsec chop 
throw was chosen in order to optimise the sky subtraction. 
The chop direction is fixed in celestial co-ordinates, lying di- 
rectly north-south in the Lockman Hole East, and 48° east 
of north in ELAIS N2, parallel to the long axis of the survey 
area (to ensure the chop remains within the boundaries). 

The preliminary observing strategy was to build up a 
linear strip in each field by overlapping SCUBA jiggle-maps, 
staggered by half the width of the array. The data collected 
from the pilot Lockman Hole East run in March 1998 were 
composed of a series of jiggle-maps, forming a deep strip at 
the centre of the survey area with crgso — 1.6 mjy/beam. In 
subsequent observing runs, the strategy was modified to that 
of an overlapping 'tripod' positioning scheme, with three 
offset grids of hexagonally-close-packed pointings, and the 
overall integration time reduced by a factor of ~ 2 to make 
more efficient use of the allocated observing time. Each tri- 
pod position received 3.0 hours of integration time, yielding 
a typical noise level of crgso — 2.5 mJy/beam. In order to give 
near homogeneous sky noise across the final coadded map, 
high and low airmasses were evenly distributed throughout. 
We have mapped a total area of 260 arcmin^, up to and 
including observations made in March 2001. 

Prior to 15th December 1999, the observations were 
made using the narrow-band 850/450 /im filters. The ma- 
jority of the observations carried out after this date made 
use of the more sensitive wide-band 850/450 /itm filters. Since 
the central wavelength of the 850 /im filter did not change by 
more than 1 /im, the data taken with the wide- and narrow- 
band filters were directly combined in producing the final 
images. 

Sky dips and pointings were carried out, on average, ev- 
ery 1.5 — 2 hrs. The observations were carried out in dry 
weather conditions, under the constraint T225GHz < 0.08. 
Primary/secondary calibration observations were taken at 
the start and end of every shift, with additional beam maps 
of the pointing source (0923-1-392 in the case of the Lockman 
Hole East, and 3C345 in the case of ELAIS N2) being made 
every 3-4 hours to provide tertiary fiux-density calibration. 
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ELAIS N2 


ELAIS N2 


Lockman Hole East 


Lockman Hole East 








Uniform noise 


Full area 


Uniform noise 


Full area 


Area (sq. arcmin) 






63 


102 


55 


122 


% Diflferential 


5 


rnjy 


12.3 ±4.6 


8.4 ±3.2 


37.1 ±10.3 


20.5 ±5.1 


completeness 


8 


rnjy 


53.8 ±10.2 


34.1 ±6.1 


68.3 ±12.9 


45.5 ±7.7 


at flux density 


11 


rnJy 


84.5 ±12.1 


62.4 ±8.2 


82.9 ± 15.4 


47.8 ± 7.3 


limit 


14 


rnjy 


98.3 ±12.8 


74.0 ±8.8 


97.7 ±15.1 


73.0 ±9.1 


% Integral 


5 


rnjy 


45.4 ±2.1 


30.8 ±1.3 


50.2 ±2.5 


31.6 ±1.4 


completeness 


8 


rnjy 


79.5 ±3.2 


57.8 ±2.1 


80.7 ±3.8 


52.2 ±2.0 


to flux density 


11 


mjy 


92.0 ±4.3 


68.8 ±2.9 


92.3 ±5.2 


62.2 ±2.8 


limit 


14 


mjy 


98.2 ±7.7 


73.0 ±5.1 


96.0 ±8.8 


72.4 ±5.2 


Boosting 


5 


rnjy 


1.28 ±0.05 


1.28 ±0.05 


1.35 ±0.06 


1.59 ±0.14 


factor 


8 


mJy 


1.16 ±0.04 


1.20 ±0.04 


1.12 ±0.04 


1.15 ±0.04 


at flux density 


11 


mJy 


1.06 ±0.03 


1.08 ±0.03 


1.10 ±0.03 


1.14 ±0.03 


limit 


14 


rnjy 


1.02 ±0.02 


1.04 ±0.02 


1.03 ±0.03 


1.05 ±0.02 


Positional 


5 


mJy 


2.42 ±0.57 


2.42 ±0.57 


4.52 ±0.53 


4.85 ± 0.49 


uncertainty 


8 


rnjy 


3.63 ±0.34 


3.92 ±0.38 


3.87 ±0.47 


3.82 ±0.40 


at flux density 


11 


mJy 


2.53 ±0.27 


2.60 ±0.23 


2.81 ±0.28 


2.96 ±0.24 


limit /arcsec 


14 


mJy 


2.39 ±0.18 


2.65 ±0.18 


2.17 ±0.17 


2.45 ±0.18 



Table 1. Results of the completeness simulations at the 3.50(t significance level, in which individual sources of the specified flux density 
given in column 2 were added into the unconvolved survey data taken up to and including 16th August 2000, and retrieved by means of 
the source extraction algorithm outlined in Section 4. 



3.5 - 




2.5 - 



nput flux /mJy 



Figure 1. The dependence of positional error on input source flux density in the range 5 — 15 mjy, based on a single template source 
repeatedly being added into and then retrieved (> 3.50(t) from the real unconvolved signal maps. The error bars give the standard error 
on the mean. Plot (a) gives the results from the uniform noise region, whereas plot (b) gives the results for the non-uniform edge regions 
which have not received the full integration time. The solid lines are a best-flt to the plotted data points. 
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Figure 2. The percentage of sources returned against input source flux density in the range 5 — 15 mJy, based on a single template 
source repeatedly being added into and retrieved (> 3.50a-) from the real unconvolved signal maps. The error bars give the standard 
error on the mean. Plot (a) gives the results from the uniform noise region, whereas plot (b) gives the results for the non-uniform edge 
regions which have not received the full integration time. The solid curves show a best-fit exponential-based model to the plotted data 
points. 



3 DATA REDUCTION AND MAP 
PRODUCTION 

The data have been reduced using two independent software 
packages; SURF (Jenness et al. 1997) and an IDL based re- 
duction package (Serjeant et al. 2001). The two mechanisms 
follow a very similar core reduction process, the main differ- 
ence being in the production of the final maps. In addition 
to the output signal maps, the IDL routines automatically 
create corresponding uncorrelated noise maps by consider- 
ing the signal variance. The final images produced by the 
SURF and IDL reductions were found to be consistent to 
within a few per cent of each other. However, given the im- 
proved knowledge of the noise local to each source, the IDL 
mechanism was adopted as the primary reduction method 
and it is this process which we describe here. 

After combining the positive and negative beams, the 
data were flat-fielded and then corrected for atmospheric ex- 
tinction. On the nights when the CSO 225 GHz opacity me- 
ter made sufficiently good measurements around the time 
of a map, each jiggle pointing of the data was extinction- 
corrected by means of a polynomial fit to the 225 GHz opac- 
ity readings, interpolated to 850 and 450 /im by means of the 
relations in Archibald et al. (2000b). If the CSO 225 GHz 
opacity meter was broken, or if the 225 GHz opacity mea- 
surements around the time of a jiggle map were poor, the 
extinction correction was assessed instead by a linear inter- 
polation between consecutive skydips. 

An iterative procedure was employed to deglitch the 



data-sets and subtract any residual sky emission not re- 
moved by chopping and nodding. Each iteration made a 
time-dependent estimate of the noise for each bolometer, re- 
moving any spikes by means of a 3cr-clip. A temporal modal 
sky level was determined by a fit to all bolometers in the 
array and subtracted from the data. With each consecutive 
iteration, the deglitching process makes a harder cut. Noisy 
bolometers were assigned a low inverse variance weight in 
this way. 

To construct the final images the signal data were 
binned into 1 arcsec pixels, creating 'zero-footprint' maps 
with a corresponding noise value determined from the sig- 
nal variance. The term 'zero-footprint' is an analogy with 
the drizzling algorithm (Fruchter & Hook, 2002). A stan- 
dard shift-and-add technique takes the flux in a given de- 
tector pixel and places its flux into the final map, over an 
area equivalent to one detector pixel projected on the sky. 
Drizzling on the other hand takes the flux and and places it 
into a smaller area in the final map. Simulations have shown 
that this helps preserve information on small angular scales, 
provided that there are enough observations to fill in the re- 
sulting gaps. The area in the coadded map receiving the flux 
from one detector pixel is termed the footprint. Our method 
is an extreme example of drizzling; we take the data from 
each 14.5 arcsecond bolometer and put the flux into a very 
small footprint (the 'zero-footprint'), one arcsecond square. 
Unlike in the standard SURF reduction, there is no intrin- 
sic smoothing or interpolation between neighbouring pixels 
in this rebin procedure. Although there is some degree of 
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(b; Non — uniform noise region 



10 
hput flux /mJy 



15 5 



10 15 

hpuf flux /mJy 



Figure 3. The factor by which the input source flux is boosted (when retrieved from the sub-mm map) plotted against input source flux 
density in the range 5 — 15 mJy, based on a single template source repeatedly being added into and retrieved (> 3.50it) from the real 
unconvolved signal maps. The error bars give the standard error on the mean. Plot (a) gives the results from the uniform noise region, 
whereas plot (b) gives the results for the non-uniform edge regions which have not received the full integration time. The solid curves 
show a best-fit exponential-based model to the plotted data points. 



correlation between pixels in the output zero-footprint sig- 
nal maps in terms of the beam pattern, the corresponding 
pixel noise values represent individual measurements of the 
temporally varying sky noise averaged over the dataset in- 
tegration time, at a specific point on the sky, and are hence 
statistically independent from their neighbours. In essence 
we have a very oversampled image with statistically inde- 
pendent pixels. Statistical non-independence would refer to 
pixel-to-pixel crosstalk, which is not the case here. A final 
4(j-clip on the signal-to-noise was carried out to remove any 
remaining 'hot pixels'. A noise- weighted convolution with a 
beam-sized Gaussian point spread function (PSF) produces 
realistic smoothed maps of the survey areas and is able to 
account for variable signal-to-noise between individual pix- 
els. 

Due to the overlap of individual jiggle-maps, and the 
fact that the observations were spread over 41 observing 
shifts, data missing due to bad bolometers are generally 
widely distributed, and only one region of the final 450 /xm 
ELAIS N2 map has been noticeably affected. This is in the 
region of RA 16:36:48.8 DEC -1-41:02:44, but is in fact a 
region which contains no significant 850 /im sources. Noise 
arising from undersampling by bolometers at the edge of 
the array has been evened out as far as possible by means 
of the tripod positioning scheme and only significantly af- 
fects the perimeter of the final mosaiced maps. Once these 
edges have been trimmed off, the resulting images have al- 
most uniform noise across the central regions of the survey 
areas (o-gso = 2.5 ±0.7 mjy/beam). The nature of the tripod 



positioning scheme does, however, lead to a noisier border 
of thickness ~ 1/2 an array width surrounding the uniform 
noise region due to incomplete observations of one or more 
overlap positions, and hence less total integration time. 



3.1 Calibration 

Dunne & Bales (2001) have examined the issue of calibra- 
tion errors in detail using the data from the 'SCUBA Lo- 
cal Universe Galaxy Survey' taken with the narrow band 
850/450 /im filters. They find that an aperture calibration 
factor (ACF) is more accurate in determining source flux 
densities than the method of determining the flux conversion 
factor (FCF) by fitting to the calibrator peaks. At 850 /im 
the standard deviation on the ACF is ~ 6%, compared with 
~ 8% when applying an FCF fitted to the peak. The im- 
provement is particularly significant at 450 /im where an 
aperture calibration method yields a typical standard de- 
viation of ~ 10% compared with ~ 18% applying a flux 
conversion factor fitted to the peak. The greater variation 
in gain when fitting to the calibrator peaks arises because of 
the sensitivity of the peak flux density to changes in beam 
shape, due to sky noise, pointing drifts, chop throw, and the 
dish shape (particularly important at 450 /im). 

Our IDL pipeline does, however, have the considerable 
advantage of producing uncorrelated noise maps, which have 
allowed us to apply a maximum-likelihood method to mea- 
sure the statistical significance of each peak in our images, 
(see Section 4, and Appendix A). Our source-extraction al- 
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(a) Uniform noise region 



(b) Non — uniform noise region 
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Figure 4. Based on the fitted flux-density boosting curves shown in the preceding Figure, these plots show the predicted output flux 
density against input flux density in the range 5 — 15 mJy. Plot (a) gives the results from the uniform noise region, whereas plot (b) gives 
the results for the non-uniform edge regions which have not received the full integration time. The dashed line in both cases delineates 
the ideal case of output flux density equal to input flux density. 



gorthim provides a simultaneous fit to all peaks in the im- 
age, and is therefore able to decouple any partially-confused 
sources whilst still recovering the additional signal-to-noise 
yielded by the negative sidelobes. The resulting best-fit 
model therefore allows us to measure the flux densities of 
our significant detections (> S.OOa) directly, as well as yield- 
ing properly-quantified errors on this value. At 850 /im the 
accuracy gained by this approach more than compensates 
for the extra ~ 2% error in the calibration and so we have 
chosen to employ a calibration system based on fitting to 
the peak values. For consistency, we have adopted the same 
approach when considering the 450 /xm data in this paper. 
The short-wavelength data will, however, be considered in 
more detail by Fox et al. (2002). 

In order to calibrate the data, each shift was divided 
in two and all available information was utilised to calculate 
the mean flux-conversion factors applicable to each half shift 
by fitting to the calibrator peaks. Uranus and Mars were ob- 
served as primary calibrators when available, with the plane- 
tary flux densities taken from the JCMT FLUXES program. 
Additionally, the secondary calibrators CRL618, OH231.8, 
CRL2688 and IRC10216 were used. Since IRC10216 is vari- 
able with a period of 635 days, it was assumed to have con- 
stant flux density over an 8/9 night observing run, calculated 
using all other calibration data for that run. Beam maps of 
the pointing sources were used for tertiary calibration. A 30- 
arcsec chop-throw was used for observing both the survey 
areas and the calibrators, in order to cancel out the effects 
of beam distortion as far as possible. Taking into account 
the errors on the calibrator flux densities, the variation in 



the flux-conversion factor over half a shift, and errors in the 
measured atmospheric opacity, the average calibration error 
over the survey maps was found to be 9% at 850 /im, and 
20% at 450 ^m. 



3.2 Pointing Problems 

During the course of the survey, 3 different pointing prob- 
lems are now known to have afflicted the JCMT. 

First, prior to 28th July 1998, a bug in the chopping 
software of SCUBA meant that the chop direction was not 
updated once an observation had begun, and consequently 
the negative sidelobes accompanying a source in this data set 
will appear smeared and therefore shallower than expected. 
However, the effect is small when compared to the sky noise 
of ~ 2.5 mjy/beam and affects only a small portion of our 
data (51/502 maps). 

Second, some of the jiggle-maps taken since June 1999 
may have been affected by an elevation drive pointing glitch 
(Coulson 2000). Observations appear to suffer a 4 arcsec 
jump in the elevation when passing through transit. Addi- 
tionally, pointing on the opposite side of the meridian to 
the survey fleld is likely to introduce a similar effect, of 
order 1.5-2 arcsec in magnitude. The fact that each area 
of sky is covered by 12 individual jiggle-maps significantly 
reduces the likelihood of a large increase in positional un- 
certainty and source smearing. Under the most pessimistic 
assumptions, a maximum of 83/502 maps could be affected, 
of which a maximum of 3/12 cover any one area of sky. 

Finally, data taken after 15th July 1999 were affected by 
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a non-synchronization of the GPS and SCUBA data acquisi- 
tion computer clocks. This problem has now been corrected 
for in the data reduction (Jenness 2000). 



4 SOURCE EXTRACTION TECHNIQUE 

The chopping-nodding mechanism of the telescope provides 
a valuable method of discriminating between real detections 
and spurious noise spikes in the data. The 30-arcsec chop 
throw is small enough to fall onto the SCUBA array and 
so negative sidelobes, half the depth of the peak flux, are 
produced on either side of a real source, 30 arcsec away 
from the position of the central peak. This side-lobe signal 
can be recovered to boost the overall signal-to-noise ratio of 
a detection. 

For well-separated sources, convolving the images with 
the beam is formally the best method of source extrac- 
tion (Eales et al. 1999, 2000, Serjeant et al. 2001). How- 
ever, following a careful examination of the reduced data it 
became clear that some of the sources were partially con- 
fused, particularly in the map of ELAIS N2 where the nega- 
tive sidelobes of individual sources had overlapped and were 
therefore somewhat deepened relative to both source peaks. 
Consequently, in order to decouple any confused sidelobes a 
source-extraction algorithm was devised based on a simul- 
taneous maximum-likelihood fit to the flux densities of all 
potentially significant peaks in the beam-convolved maps. 

Using a peak-normalised beam map as a source tem- 
plate (CRL618 in the case of the Lockman Hole East, or 
Uranus in the case of ELAIS N2), a basic model was con- 
structed by centring a beam map at the positions of all peaks 
> 3 mjy as found in the Gaussian-convolved maps. The 
heights of each of the positioned beam maps were then cal- 
culated such that the final multi-source model provided the 
best description of the sub-mm sky, as judged by a minimum 
X'^ calculation made feasible by the independent data-points 
and errors yielded by the zero-footprint IDL-reduced maps. 
A full mathematical description of this technique is provided 
in Appendix A. 

The actual maps and source analysis are presented in 
Section 6. First, however, it is necessary to understand the 
selection effects that arise in an analysis of this sort. 



5 SIMULATIONS 

In order to assess the effects of confusion and noise on 
the reliability of our source-extraction algorithm, Monte 
Carlo simulations were carried out using the 220 out of 
260 arcmin^ of data taken by the end of August 2000. The 
dependences of positional error, completeness, and error in 
reclaimed flux density, on source flux density and noise in 
the maps were determined by planting individual sources of 
known flux density into the real ELAIS N2 and Lockman 
Hole East SCUBA maps. This has the advantage of allow- 
ing us to test the source-reclamation process against the real 
noise properties of the images. However, these simulations 
do not allow assessment of the level to which false or con- 
fused sources can contaminate an extracted source list. We 
have therefore also created a number of fully-simulated im- 
ages of the survey areas by assuming a reasonable 850 pun 



source-counts model. The results of analyzing these two sets 
of simulations are discussed in the following two subsections. 

5.1 Simulations building on the real survey data 

A normalised beam map was used as a source template. At 
0.5 mJy flux intervals, from 5 to 15 mJy, individual sources 
were added into the unconvolved zero-footprint signal maps. 
This was done one fake source at a time, so as not to en- 
hance significantly any existing real confusion noise within 
the image. The source-extraction algorithm was then re- 
run. This exercise was repeated for 100 different randomly- 
selected positions on each image, at each flux level, so that 
source reclamation could be monitored as a function of input 
flux and position/noise-level within the maps. Source recla- 
mation was deemed to have been successful if the source- 
extraction algorithm returned a source with signal-to-noise 
> 3.50 within less than half a beam width of the input posi- 
tion. When analysing the output data, ~ 10% of the input 
sources were found to lie within half a beam width of an- 
other peak > 5 mJy already present in the image. Any such 
source was excluded from statistical analyses of positional 
error, completeness and error in reclaimed flux density, since 
the proportion of blends produced by two sources, both with 
comparably bright fluxes in the region of 8 mJy, should be 
negligible within our sample. 

Figs 1 to 4 show the dependences of mean positional 
error, percentage of sources returned, systematic offset in 
returned flux density, and predicted output against input 
ffux density, for both regions of uniform and non-uniform 
noise. The error bars show the standard error on the mean. 
The solid lines in Fig 1, and solid curves in Figs 2 and 3, 
are best-fit models to the plotted data points. From Fig 3, 
we see that the effect of noise and confusion is to produce 
systematic 'flux boosting', the retrieved flux density always 
being greater than the input value. The solid curves in Fig 
4 show the predicted output versus input flux density based 
on the fitted flux-density boosting models of Fig 3. We find 
that there is a strong correlation with input flux density for 
each of these quantities, the scatter being somewhat larger 
in the areas of non-uniform noise lying towards the edge, 
as would be expected. The positional error and flux-density 
boosting factors decrease with increasing input flux density, 
whilst the percentage of sources returned at > 3.50(t in- 
creases. The effect of additional noise in the non-uniform 
regions considerably increases the uncertainty in positional 
error, as well as the mean flux-density boosting factor. It also 
serves to typically half the percentage of sources reclaimed. 
A summary of the results at 5, 8, 11 and 14 mJy is given 
in Table 1. The integral completeness down to each flux- 
density limit has been calculated by a counts-slope- weighted 
sum of the fraction of sources returned at each specific flux- 
density level, based on an adopted 850 /im source counts 
model consistent with the measured source counts, given in 
Section 5.2. As detailed in Table 1, the typical complete- 
ness at 8 mJy is ~ 80% within the uniform noise regions of 
the images, with flux-density boosting factors of ~ 15%. By 
11 mJy completeness has risen to ~ 90%, while flux-density 
boosting has dropped to < 10%. 
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Figure 5. A pure noise image of the ELAIS N2 region, created from the real data by corrupting the astrometry of the individual jiggle 
maps to produce a 'scrambled' map. This image includes the jiggle maps taken up to and including 16th August 2000. There are no 
source detections at > 3.50cr. 
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Figure 6. Fully-simulated ELAIS N2 SCUBA 850 /im map, using the adopted counts model given in Section 5.2. The image noise was 
simulated using the scrambled map shown in Figure 5. 
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Figure 7. A pure noise image of the Lockinan Hole East region, created from the real data by corrupting the astrometry of the individual 
jiggle maps to produce a 'scrambled' map. This image includes the jiggle maps taken up to and including 16th August 2000. There is 
one spurious source detection at > 3.50cr which is circled on the image. 
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Figure 8. Fully-simulated Lockman Hole East SCUBA 850 /im map, using the adopted counts model given in Section 5.2. The image 
noise was simulated using the scrambled map shown in Figure 7. 
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78.6 ±16.8 


63.0 ±11.7 
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5.6 
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% count correction factor 


88.9 ±22.2 


112.0 ±21.2 


75.0 ±21.7 


90.0 ±21.2 


82.4 ±15.6 


102.2 ±15.1 



Table 2. Results of the completeness simulations at the 4.00(T significance level, in which each field was fully simulated 4 times using an 
adopted counts model consistent with the measured source counts. The sources were then retrieved by means of the maximum-likelihood 
estimator outlined in Section 4. 
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% completeness 


81.3 ±22.5 


60.7 ±14.7 
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77.8 ± 20.8 


85.7 ±17.5 


67.4 ±12.1 


% confused/random 


17.4 ± 
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15.6 ± 7.0 


21.1 ±10.5 


20.8 ± 9.3 


19.0 ± 
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17.9 ± 5.6 


% count correction factor 


69.6 ±17.4 


87.5 ±16.5 


63.0 ±18.2 


75.0 ±17.7 


66.7 ±12.6 


82.1 ±12.1 



Table 3. Results of the completeness simulations at the 3.50(T significance level, in which each field was fully simulated 4 times using an 
adopted counts model consistent with the measured source counts. The sources were then retrieved by means of the maximum-likelihood 
estimator outlined in Section 4. 
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% completeness 


81.3 ±22.5 


67.9 ±15.6 


91.7 ±27.6 


77.8 ±20.8 


85.7 ±17.5 


71.7 ±12.5 


% confused/random 


26.9 ±10.2 


40.0 ± 8.9 


51.6 ±12.9 


57.1 ± 10.8 


40.4 ± 


8.4 


48.5 ± 7.0 


% count correction factor 


61.5 ±15.4 


56.0 ±10.6 


38.7 ±11.2 
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Table 4. Results of the completeness simulations at the 3.00(T significance level, in which each field was fully simulated 4 times using an 
adopted counts model consistent with the measured source counts. The sources were then retrieved by means of the maximum-likelihood 
estimator outlined in Section 4. 
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5.2 Completely simulated maps 

Adopting an 850 /im source counts model consistent with 
the measured point source counts, and which does not over- 
predict the FIR background, we have created fully-simulated 
versions of our two survey fields. In order to produce a real- 
istic model of the background counts, the number of sources 
expected at 0.1 mjy intervals was Gaussian randomised and 
then the appropriate sources were planted into the simu- 
lated data sets at random positions. The whole image was 
then convolved with the beam. At all flux-density levels the 
sources were assumed to be un-clustered (although see Sec- 
tion 9.4 and Almaini et al. 2002). 

We then constructed a pure noise image for each of our 
two fields by scrambling the astrometry of the individual 
SCUBA observations used to produce the maps. These may 
be seen in Figs 5 and 7. After creating the scrambled noise 
maps we first ran the source-extraction algorithm on these to 
assess the frequency with which spurious sources would ar- 
tificially appear on the basis of random statistics. The noise 
images yielded only a single artificial source at > 3.50a, 
and a total of 12 spurious sources at > 3.00a, which is 
not unexpected given the ~ 1200 beams across the areas 
of these maps. Adding the simulated background counts to 
the scrambled maps produces the final simulated images. 
The source-extraction algorithm was applied as previously. 

Four realisations of the sky in each field were generated 
for the source counts model 



dN{> s) ^ m f(S_Y/^\>' 

dS So \\So) \So 



(1) 



where iV,) = l.lxlO^deg^Vjy"^ , S'o = 1.79 mJy, a = 0.0, 
and P — 3.13. Examples of the simulated ELAIS N2 and 
Lockman Hole East maps for our adopted counts model may 
be seen in Figs 6 and 8. 

Tables 2 to 4 summarise the results at 4.00, 3.50 and 
3.00a respectively for the simulated fields at 8 mJy. In these 
tables and the discussion that follows, we have adopted the 
practical criterion that any source retrieved at the 8 mJy 
level which ultimately could not be identified with a sin- 
gle source brighter than 5 mJy (and therefore could not be 
confirmed through interferometric follow-up for example), 
should be regarded as spurious or the result of confusion. 
In the uniform noise regions, contamination from spurious 
sources or confused sources < 5 mJy is ~ 40% at S/N > 3.00, 
dropping to ~ 20% at S/N > 3.50 and < 10% at S/N > 4.00. 
The completeness at a measured flux density of 8 mJy is 
~ 85% within the regions of uniform noise, consistent with 
the completeness determined from the simulations using the 
real survey data. We also flnd that although the real source 
counts (ie. total number of sources out minus the number 
of spurious sources) at > 3.50a agree well with the num- 
ber planted into the simulations at 8 mJy, this is partly by 
chance; the number of sources boosted from a lower flux- 
density level, between 5 and 8 mJy, is approximately the 
same as the number of real sources > 8 mJy which are 
not retrieved by the source-extraction algorithm. By low- 
ering our signal-to-noise threshold from 3.50 to 3.00, there 
is very little improvement (if any) in completeness, however 
the proportion of spurious sources detected is much more 
serious. These results provide a clear motivation for setting 
a signal-to-noise threshold of > 3.50 when determining the 



source number counts as opposed to the > 3.00a level ap- 
plied in some other surveys. Our results are consistent with 
those of Hughes & Gaztanaga (2000), who have simulated 
sub-millimetre galaxy surveys for a range of wavelengths, 
spatial resolutions and flux densities. 



6 SOURCES 

Sources down to a significance level of > 3.00a at 850 fim are 
listed in Table 5 for ELAIS N2 and Table 6 for the Lockman 
Hole East, ranked in order of formal significance as deter- 
mined from our maximum-likelihood extraction technique. 
The top section of each table, containing 7 and 12 sources in 
ELAIS N2 and Lockman Hole East respectively, is confined 
to sources for which the S/N > 4.00, and the middle section 
to sources for which S/N > 3.50. Since our simulations in- 
dicate that 70-80% of these sources are expected to be real 
individual sources with flux densities within ~ 25% of the 
extracted flux density, we recommend that follow-up obser- 
vations are confined to these abbreviated S/N > 3.50 source 
lists. The 850 /im flux densities in column 4 are the best- 
fit scaling factors determined by our extraction algorithm, 
with the quoted errors also incorporating the 9% calibration 
error. 

Figs 9 and 10 show the beam-sized Gaussian convolved 
survey maps at 850 /im for ELAIS N2 and the Lockman 
Hole East respectively, with the sources > 3.50a circled and 
numbered. We have also applied our maximum-likelihood 
estimator to the 450 /im data (Fox et al. 2002), centring 
normalised 450 /im beam maps at the same positions as in 
the 850 /im images. Column 7 gives the 450 fim flux densities 
if there is a detection in the short- wavelength data with S/N 
> 3.00a, or, more commonly, a 3a- upper-limit. We also car- 
ried out a search for peaks in the short- wavelength data lying 
within a search radius of 8 arcsec from the 850 /im source 
peaks (the combined result of adding in quadrature the 1/2 
beam- width positional uncertainties of the long- and short- 
wavelength data) and again employed our source-extraction 
algorithm. The flux densities of any detections with S/N 
better than 3.00a are given in column 9, with their signifl- 
cances and positional offset from the long-wavelength posi- 
tion in columns 10 and 11. The errors quoted on the 450 /im 
flux densities also include the 20% calibration error. 

Given that the FWHM of the JCMT beam at 850 ^m is 
14.5 arcsec, the source positions are only likely to be accu- 
rate to within ~ 3 — 4 arcsec. Consistent with this, our sim- 
ulations suggest that an 8 mJy source, detected at > 3.50a, 
should have accurate astrometry to ~ 3.5 arcsec in uniform 
noise regions and ~ 4 arcsec in the non-uniform noise areas, 
the positional uncertainty decreasing with increasing flux. 
However, these simulations do not account for the elevation 
drive pointing glitch mentioned in Section 3.2, and so it is 
possible that the astrometric accuracy of a source affected 
by this problem could be somewhat worse than this. 

Analysis of our fully-simulated survey areas suggests 
that of the 19 sources with S/N > 4.00, only one is likely 
to be spurious. At the 3.50a threshold we have 38 sources, 
of which ~ 7 may be spurious. The situation is markedly 
worse at 3.00a, where ~ 30 of the supposedly-significant 
peaks may in fact be due to noise or confusion of low flux- 
density sources. 
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Figure 9. The 850 ^m image of tlie ELAIS N2 field, smootlied witli a beam-size Gaussian (14.5 arcsec FWHM). The numbered circles 
highlight those sources found at a significance of > 3.50, the labels corresponding to the numbers in Table 5. 



In the short-wavelength data, we can find only 3 detec- 
tions at > 4.00(7 within an 8 arcsec radius of the 850 /im 
positions. These 450 /im detections appear to be associ- 
ated with the most significant Lockman Hole East 850 /im 
source and the two most significant ELAIS N2 850 /im 
sources. There are additionally another five 450 ^m detec- 
tions at better than S.OOcr, of which all but one is associ- 
ated with a > 3.50(7 850 /xm source. We must be cautious 



when considering the credibility of the lower significance 
(> 3.00(7) 450 jj,m detections. The atmospheric opacity at 
450 ^m is much more sensitive to changes in the amount 
of water vapour present than at 850 /im, and so the re- 
sulting short-wavelength maps are noisier than their long- 
wavelength counterparts. Furthermore, the smaller beam- 
size (7.5 arcsec FWHM at 450 /xm as opposed to 14.5 arcsec 
FWHM at 850 /im) means that the short-wavelength maps 
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12.9 ±3.9 


3.51 


non-uni 


< 18 










17 


16:36:51.4 


41:05:06 


5.7 ±1.7 


3.50 


uniform 


< 23 










18 


16:36:11.4 


40:59:26 


20.8 ±6.2 


3.49 


non-uni 


< 138 










19 


16:36:35.9 


41:01:38 


9.4 ±2.8 


3.49 


uniform 


< 76 










20 


16:36:49.4 


41:04:17 


6.6 ±2.0 


3.48 


uniform 


22 ±8 


3.32 


20 ± 8 


3.05 


3.2 


21 


16:37:10.5 


41:00:49 


4.9 ±1.5 


3.47 


uniform 


<23 










22 


16:36:27.9 


40:54:04 


13.4 ±4.1 


3.46 


non-uni 


< 102 










23 


16:37:19.5 


41:01:38 


12.4 ±3.8 


3.45 


non-uni 


<40 










24 


16:36:26.9 


41:02:23 


10.4 ±3.2 


3.43 


non-uni 


< 44 










25 


16:36:18.3 


40:59:12 


12.1 ±3.8 


3.37 


non-uni 


<57 










26 


16:36:32.0 


41:00:05 


9.7 ±3.1 


3.26 


uniform 


< 27 










27 


16:36:48.3 


41:03:52 


6.6 ±2.1 


3.25 


uniform 


< 21 










28 


16:36:47.2 


41:04:48 


6.3 ±2.0 


3.24 


uniform 


< 22 










29 


16:36:24.1 


40:59:35 


9.9 ±3.3 


3.14 


uniform 


<70 










30 


16:37:07.5 


41:02:37 


4.6 ±1.5 


3.13 


uniform 


<40 










31 


16:36:28.1 


41:01:41 


6.8 ±2.3 


3.07 


uniform 


<36 










32 


16:36:39.8 


41:00:34 


5.4 ±1.8 


3.07 


uniform 


<27 










33 


16:36:50.5 


40:58:54 


4.8 ±1.6 


3.06 


uniform 


<34 










34 


16:36:27.0 


40:58:15 


6.9 ±2.3 


3.05 


uniform 


<36 










35 


16:36:44.8 


40:56:51 


5.5 ±1.9 


3.02 


uniform 


<23 










36 


16:36:52.9 


41:02:52 


3.8 ±1.3 


3.00 


uniform 


n/a 











Table 5. Sources retrieved in the long-wavelength ELAIS N2 image above the 3.00iT significance level using the source-extraction 
algorithm outlined in Section 4. We used a peak-normalised calibrator beam-map as a source template, and constructed a basic model 
by centring a beam-map at the positions of all peaks > 3 mJy as found in the 14.5" FWHM Gaussian convolved maps. The heights of 
each of the positioned beam maps were allowed to vary simultaneously such that the final model satisfied a minimised x^ fit. We suggest 
that follow-up observations be confined to those detections > 3.50(7 given the serious problem of contamination by spurious sources at 
the 3.00(7 level. Calibration errors of 9% and 20% at 850 /im and 450 /im respectively, are included in the error on the flux densities. 



will contain a greater number of spurious detections. For 
S/N > 3.00 we would expect 53 spurious sources across the 
full 260 arcmin^ of short-wavelength data ( and 14 spurious 
detections in the long-wavelength data) based on Gaussian 
statistics alone. 



A detailed discussion of individual sources may be found 
in paper II of the SCUBA 8 mJy survey (Fox et al., 2002), 
along with a multiwavelength analysis considering possi- 
ble identifications and SED-based redshift constraints. The 
450 /jm images of the survey fields are also presented by Fox 
et al. (2002). 



7 NUMBER COUNTS 

In order to obtain the best estimate of the number counts 
at 850 /im, we must account for the effects of flux-density 
boosting, incompleteness and contamination from spurious 
sources. Using the results of our simulations carried out 
on the real survey data, the observed flux densities of the 
> 3.50iT sources were corrected statistically using the out- 
put flux density versus input flux density plots shown in 
Fig 4, for the appropriate image noise region. The number 
of boost-corrected sources in bins 0.5 mJy wide, centred on 
5.0, 5.5, 6.0 mJy and so on up to 15 mJy, were then amended 
for completeness using the best fit results of the real-data 
simulations shown in Fig 2. The integrated source counts 
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Figure 10. The 850 ^m image of ttie Lockinan Hole East field, smootlied witti a beam-size Gaussian (14.5 arcsec FWHM). The numbered 
circles highlight those sources found at a significance of > 3.50, the labels corresponding to the numbers in Tabic 6. 



were obtained by summing the number of sources in the bins 
down to the required threshold flux. Since the scatter in the 
simulation results is substantially higher in the non-uniform 
noise regions (reflected by reduced x^ values of the plotted 
best-fit curves being a factor of ~ 5 higher in the edge re- 
gions), we have based our estimation of the 850 /um source 
counts purely on the sources detected in the areas of uni- 
form noise, covering an area of 190 arcmin'^. In practice this 



means that 5/38 of the sources above the 3.50(t significance 
threshold have been excluded. This process is summarised 
in Table 7, with the raw and amended counts per square 
degree shown in Table 8, along with Poisson errors. 

The applied correction factors have been determined 
from simulations carried out on the real images and so are in- 
dependent of source-count models. We cannot, however, at- 
tempt to correct for the presence of random detections with- 
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Source 


RA 


DEC 


•5850 


S/N 


Noise 


541 


50 or 3(7 


S/N 5450 of any 


S/N 


Distance 




{J2000) 


{J2000) 


/mjy 


850 ^lm 


Region 


up 


. lim. at 


450 /im > 3(7 source 


450 ^lva 


from 














850 /im pos. 


at < 8.0" from 


at 


850 Atm 














/mJy 


850 /im 850 /im pos 


450 ura 


peak 


















pos /mJy 


pos 


/arcsec 


01 


10:52:01.4 


57:24:43 


10.5 ±1.6 


8.10 


deep 


< 


43 


45 ±13 


4.64 


5.0 


02 


10:52:38.2 


57:24:36 


10.9 ±2.4 


5.22 


uniform 


< 


41 








03 


10:51:58.3 


57:18:01 


7.7 ±1.7 


5.06 


deep 


< 


54 


39 ±14 


3.38 


3.2 


04 


10:52:04.1 


57:25:28 


8.3 ±1.8 


5.03 


deep 


< 


51 








05 


10:51:59.3 


57:17:18 


8.6 ±2.0 


4.57 


deep 


< 


34 








06 


10:52:30.6 


57:22:12 


11.0 ±2.6 


4.50 


uniform 


< 101 








07 


10:51:51.5 


57:26:35 


8.1 ±1.9 


4.50 


deep 


< 


86 








08 


10:52:00.0 


57:24:21 


5.1 ±1.3 


4.38 


deep 


< 


17 








09 


10:52:22.7 


57:19:32 


12.6 ±3.2 


4.20 


uniform 


< 


47 








10 


10:51:42.4 


57:24:45 


12.2 ±3.1 


4.18 


uniform 


< 124 








11 


10:51:30.6 


57:20:38 


13.5 ±3.5 


4.09 


non-uni 


< 


90 


65 ±24 


3.32 


6.0 


12 


10:52:07.7 


57:19:07 


6.2 ±1.6 


4.01 


deep 


< 


55 








13 


10:51:33.6 


57:26:41 


9.8 ±2.8 


3.69 


uniform 


< 


90 








14 


10:52:04.3 


57:26:59 


9.5 ±2.8 


3.61 


uniform 


< 274 








15 


10:52:24.6 


57:21:19 


11.7 ±3.4 


3.60 


uniform 


< 


53 








16 


10:52:27.1 


57:25:16 


6.1 ±1.8 


3.56 


uniform 


< 


39 








17 


10:52:16.8 


57:19:23 


9.2 ±2.7 


3.55 


uniform 


< 


57 








18 


10:51:55.7 


57:23:12 


4.5 ±1.3 


3.55 


deep 


< 


36 








19 


10:52:29.7 


57:26:19 


5.5 ±1.6 


3.54 


uniform 


< 


39 








20 


10:52:37.7 


57:20:30 


10.3 ±3.1 


3.51 


non-uni 


< 


73 








21 


10:52:01.7 


57:19:16 


4.5 ±1.3 


3.50 


deep 


< 


35 








22 


10:52:05.7 


57:20:53 


4.7 ±1.4 


3.49 


deep 


< 


52 








23 


10:51:47.0 


57:24:51 


7.4 ±2.2 


3.48 


uniform 


< 


60 








24 


10:51:42.9 


57:24:12 


11.4 ±3.4 


3.47 


uniform 


< 


63 








25 


10:52:36.0 


57:18:20 


12.7 ±3.8 


3.46 


non-uni 


< 137 








26 


10:52:27.3 


57:19:06 


8.3 ±2.6 


3.39 


uniform 


< 


59 








27 


10:51:53.8 


57:18:47 


5.5 ±1.7 


3.38 


deep 


< 


33 








28 


10:52:34.6 


57:20:02 


10.2 ±3.2 


3.31 


uniform 


< 


74 








29 


10:52:16.4 


57:25:07 


6.7 ±2.1 


3.30 


uniform 


< 


28 








30 


10:52:42.2 


57:18:28 


11.2 ±3.6 


3.25 


non-uni 


< 


92 








31 


10:52:03.9 


57:20:07 


4.0 ±1.3 


3.24 


deep 


< 


22 








32 


10:52:00.0 


57:20:39 


4.3 ±1.4 


3.22 


deep 


< 


24 








33 


10:51:33.8 


57:19:29 


8.1 ±2.6 


3.20 


uniform 


< 


75 








34 


10:52:09.9 


57:20:40 


8.5 ±2.8 


3.16 


uniform 


< 


77 








35 


10:51:57.6 


57:26:03 


6.7 ±2.3 


3.02 


uniform 


< 


91 








36 


10:52:03.5 


57:16:54 


8.0 ±2.8 


3.00 


uniform 


< 


29 









Table 6. Sources retrieved in the Lockman Hole East field above the 3.00(T significance level using the source extraction algorithm 
outlined in Section 4. We used a peak-normalised calibrator beam-map as a source template, and constructed a basic model by centring 
a beam-map at the positions of all peaks > 3 mJy as found in the 14.5" FWHM Gaussian convolved maps. The heights of each of 
the positioned beam maps were allowed to vary simultaneously such that the final model satisfied a minimised x^ fit- We suggest that 
follow-up observations are confined to those detections > 3.50(T given the serious problem of contamination by spurious sources at the 
3.00(7 level. Calibration errors of 9% and 20% at 850 ^m and 450 ^m respectively, are included in the error on the flux densities. 



out making some assumptions about the background counts. 
Consequently, the estimated number of spurious sources is 
included in the lower error bars combined in quadrature with 
the Poisson error, rather than being applied directly to the 
actual tabulated source-count values. 

Fig 11 shows our corrected source counts, at 1 mJy in- 
tervals, from 5 to 12 mJy, along with data points from other 
SCUBA surveys (Barger, Cowie & Sanders 1999, Blain et 
al. 1999, Eales et al. 2000, Hughes et al 1998, and Small et 
al. 1997). The upper solid and dashed lines in Fig 11 are 
source-count models, derived by interpolating the 60 pm lu- 
minosity function (Saunders et al., 1990) to 850 //m, with 
an assumed dust temperature of 40K and dust emissivity in- 
dex /3 = 1.2 (consistent with measurements of local ULIRGs 



such as Arp220 eg. Dunne et al. 2000) , and assuming pure lu- 
minosity evolution of the form L{z) = L(0)(l+2:) to 2 = 2.0 
(beyond which the luminosity function is simply frozen). 
The solid line assumes an Einstein-de Sitter cosmology, 
whereas the dashed line assumes the cosmology Q.m ~ 0.3 
and ^A = 0.7 (both adopting Ho = 67 kms"^Mpc"^). The 
lower solid and dashed lines represent the same interpolated 
60 fim luminosity function with no luminosity evolution in- 
cluded, again for cosmologies Qm = 1-0, flA ~ 0.0, and 
f^M = 0.3, Oa ~ 0.7 respectively. The dotted line shows 
the adopted source counts model, which was assumed when 
creating the fully-simulated images. The dot-dash and dot- 
dot-dash lines represent a more physical form of luminos- 
ity evolution (Jameson 2000, Small et al. 2001), which is 
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Flux density 


Number 


N(>S) 


Number after 


Number after 


N{>S) 


/mjy 


observed 


observed 


flux boost 
correction 


flux boost & 

completeness 

correction 


corrected 


5.0 


3 


31 





0.00 


32.9 


5.5 


4 


28 


1 


3.59 


32.9 


6.0 


2 


24 





0.00 


29.3 


6.5 


1 


22 


3 


6.55 


29.3 


7.0 


1 


21 


3 


5.66 


22.7 


7.5 


1 


20 





0.00 


17.1 


8.0 


2 


19 


3 


4.64 


17.1 


8.5 


4 


17 


1 


1.44 


12.4 


9.0 


3 


13 


1 


1.36 


11.0 


9.5 


1 


9 





0.00 


9.6 


10.0 


1 


9 


3 


3.75 


9.6 


10.5 


2 


8 


2 


2.42 


5.8 


11.0 


3 


6 





0.00 


3.4 


11.5 


1 


3 


1 


1.15 


3.4 


12.0 


1 


2 


1 


1.13 


2.2 


12.5 


1 


1 


1 


1.11 


1.1 


13.0 











0.00 


0.0 



Table 7. A breakdown of the raw and corrected source counts at each 0.5 mJy interval bin for the sources with S/N > 3.50 from both 
survey maps, excluding those uncovered in the non-uniform noise regions. Column 1 gives the flux density, and columns 2 and 3 the 
raw number counts in the 0.5 mJy wide bin and cumulative to that flux-density level. Column 4 gives the number in each bin once 
flux-density boosting has been accounted for using the fitted curve in Fig 4(a), and similarly column 5 gives the number in each bin once 
completeness has also been taken into account using the fitted curve in Fig 2(a). Column 6 gives the cumulative corrected counts down 
to each fiux-density level. 



Flux density 


Raw 850 fivci source counts 


Corrected 850 /im source counts 


/mJy 


N{> 5)deg-2 


N{> 5)deg-2 


5 
6 
7 
8 
9 


580 ± 100 
450 ± 90 
400 ± 90 
360 ± 80 
250 ± 70 


62011J" 

ccn+lOO 
4on+90 

320t?»o 
200+^° 


10 


170 ± 60 


180l«^ 


11 


110 ± 50 


601^3^ 


12 


40 ±30 


401^0 



Table 8. The 850 fim source counts per square degree based on sources with S/N > 3.50 in both survey maps, and excluding those 
detected in the non-uniform noise regions. Column 1 gives the fiux density and column 2 the cumulative raw counts per square degree 
with the Poisson error. Column 3 gives the cumulative corrected counts per square degree, the upper error corresponding to the Poisson 
error, and the lower error accounting for both the Poisson error and the presence of spurious sources based on the simulation data. 



compatible with models of cosmic chemical evolution and 
includes a peak in the evolution function: 



L{z) = i(0)(l + 2)^''^sech^[b ln(l + z) - c]cosh^ 



(2) 



where b = 2.2 ± 0.1 and c = 1.84 ± 0.1. The dot-dash line 
assumes a cosmology flM = 1.0, f^A ~ 0.0, and the dot-dot- 
dash line Qm = 0.3, ^a = 0.7. A dust emissivity index of 
(5 = 1.2 and a dust temperature of 37K were also assumed 
(Small et al. 2001). 



8 2-POINT CORRELATION FUNCTIONS 

In order to obtain a first measure of the clustering properties 
of each of our two fields, we have used our > 3.50cr sources 
to calculate 2-point angular correlation functions for both 



ELAIS N2 and the Lockman Hole East. We created a cata- 
logue for each survey area, containing ~ 40000 fake sources 
placed according to a Poisson distribution within the bound- 
aries of the real data. Although the positions were generated 
at random, we used our uncorrelated noise maps to weight 
the number density of sources across the image, since we 
would expect to uncover a larger density of sources above 
a specified signal-to-noise threshold in areas of lower noise. 
In practise, we divided the full image into a series of sub- 
images, 20 arcsec by 20 arcsec in size, and calculated the 
mean noise level in each of these grid sections. For sources 
brighter than ~ 5 mJy, the number of sources expected 
above a constant signal-to-noise threshold will increase ap- 
proximately as A*'(> S) OC 5'"^'^, and consequently we would 
expect the relative number of sources in each sub-image to 
scale as (noise)""*"'^. 
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The angular correlation function 'w[6) describes the ex- 
cess probability of finding a source at an angular separation 
9 from a source selected at random, over a random distribu- 
tion. There are a variety of possible estimators for w{9) as 
a function of pair-count ratios. Here we have used: 



i + w{e) 



2m {DD) 
(n - 1) \dR) 



(3) 



where there are n and m objects in the real and mock cata- 
logues respectively, DD is the number of distinct data pairs 
in the real image within a bin covering a specified range of 
9, and DR is the number of cross-pairs between the real and 
fake catalogues within the same range of 9. 

The angular correlation functions may be seen with 
Poisson errors in Figs 12 and 13, where a bin size of twice 
the beam (29.0 arccconds) has been used. 



9 DISCUSSION 

9.1 Comparative Source Counts and Models 

Although our survey was undertaken with the aim of con- 
straining the brighter end of the 850 pan number counts, we 
have in fact achieved the most accurate determination to 
date of the sub-mm source counts down to Sgso — 5 mjy. 
From Fig 11 we see that at 5 mJy there is good agree- 
ment between our counts 620!tigQdeg~'^ and the values of 
Bales et al. (2000) 500 ± 200deg"^ and Barger, Cowie & 
Sanders (1999) 610tigodeg"^, to well within the quoted er- 
rors. However, at larger flux densities, the number of sources 
per square degree in both the Eales and Barger survey areas 
is somewhat lower than our counts. For example at 6 mJy 
the number count of Eales et al. is a factor of ~ 3 lower 
than our value, and that of Barger, Cowie & Sanders (1999) 
is lower by a factor of ~ 2.5. The most likely explanation 
for this discrepancy is that it is due to small number statis- 
tics rather than a real steepening of the source counts at 
> 6 mJy. We would in fact expect this to be the case if 
SCUBA sources were generally clustered on scales of a few 
arcmin (as indicated by our two-point autocorrelation func- 
tions, albeit with large errors, shown in Section 8 and dis- 
cussed in Section 9.4). The "SCUBA 8 mJy Survey" covers 
260 arcmin^ of sky, 5 times the area of the 'Canada UK 
Deep Sub-millimetre Survey (CUDSS) 14h Field' of Eales 
et al., and ~ 2.5 times the area of Barger et al., making 
us less susceptible to arcmin-scale clustering and thus leav- 
ing us much better able to constrain the number counts at 
brighter flux densities. At 8 mJy, Blain et al. (1999) suggest 
a rather higher count of 800 ± 600deg^^ than our value of 
320l]^QQdeg~^, however the large error on the Blain et al. 
number count renders their result almost meaningless. At 
flux densities > 10 mJy our survey also inevitably becomes 
limited by the lower surface density of brighter sources. In 
a similar manner to the source counts of Eales et al. and 
Barger et al. beyond 6 mJy, our counts at 11 and 12 mJy sug- 
gest a steepening of the counts slope towards brighter fluxes. 
Borys et al. (2002) report a number count of 164lggdeg~'^ 
at 12 mJy, suggesting that this is again a statistical rather 
than a real effect. 

The source counts accumulated across all the sub- 
millimetre surveys (as shown in Fig 11) are consistent with 
a population of high-redshift (z > 1), dusty, star-forming 



galaxies, analogous to the local ULIRGs, with dust temper- 
atures Td ^ 40K, emissivity /3 ~ 1.2, and pure luminosity 
evolution of the form (1 -I- z)^ out to a redshift z — 2 (and 
simply frozen beyond). We also find that our source counts 
are in very good agreement with the more physical form 
of luminosity evolution of Jameson (2000) and Small et al. 
(2002), assuming a dust temperature of 37 K, particularly 
for the cosmology Qm = 0.3, flA ~ 0.7. The lower solid 
and dashed lines show the predicted source counts using the 
same luminosity function but without any luminosity evolu- 
tion. Such a model is inconsistent with the measured number 
counts by 2 — 3 orders of magnitude, which simply serves to 
re-emphasise the extent of the evolution in the starburst 
population uncovered by SCUBA. 



9.2 Dust Masses and Star Formation Rates 

The photometric redshift estimates of the 19 sources dis- 
cussed in Fox et al. (2002) and Lutz et al. (2001) are all 
indicative of redshifts z > 1, and in the vast majority of 
cases favour z > 2. Even without knowing a precise red- 
shift, this fact coupled with an almost uniform sensitivity 
to objects between redshifts z = 1 and 2 = 8, allows us 
to make an estimate of the dust-enshrouded SFRs and dust 
masses. Assuming a dust temperature of 40K and emissivity 
index l3 = 1.2, we have interpolated the inferred 850 ^m lu- 
minosity to 60 jj,m using an optically-thin greybody approx- 
imation. The far-infrared luminosity provides a measure of 
the current SFR of massive stars. In regions of intense star 
formation, dust is heated primarily by the embedded O and 
B stars which evolve rapidly and dispense their surround- 
ing material on similarly short time-scales (~ 10^ yr, Wang 
1991), hence the rate of dust production is proportional to 
the star-formation rate. We may calculate the SFR to within 
a factor of a few, by means of 






(4) 



The value of e is uncertain to within a factor of ~ 3 (Scoville 
& Young 1983, Thronson & Telesco 1986, Rowan-Robinson 
et al. 1997, Rowan- Robinson 2000), the main sources of error 
arising from uncertainties in the initial mass function (IMF), 
the fraction of optical/UV light absorbed by the dust, and 
the time-scale of the burst. We have assumed a value of 
e = 2.1 (Thronson & Telesco 1986), which incorporates a 
burst of 0,B, and A type star-formation over ~ 2 x 10® yr, 
and assumes a Salpeter IMF. It was also assumed that all 
of the optical/UV radiation was absorbed and re-radiated 
thermally by the surrounding dust. 

In calculating the mass of dust present in each source, 
we assumed the sub-millimetre continuum to be the result of 
optically-thin thermal emission from the heated dust grains, 
with no additional contribution from bremsstrahlung or syn- 
chrotron radiation. We may then determine the dust mass 
Md by means of the relation 



Md^ 



OohsDi^ 



1 + z fcJ;'=='B(iy''«=STd) 



(5) 



where z is the redshift of the source. Sobs is the observed 
flux density, Dl is the luminosity distance, fcjj'^''* is the rest- 
frequency mass absorption coefficient, and B{y^^^^ ,Td) is 
the rest-frequency value of the Planck function from dust 
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^ Small et al. 1997 , 

o Hughes et al. 1998 \ 
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* Borys et al. 2002 
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Flux/mJy 



Figure 11. Cumulative 850 ^m source counts. The solid diamonds show the new results from the 'SCUBA 8 mjy Survey' derived in this 
paper. The upper solid and dashed curves are predicted number counts, based on the 60 ^m luminosity function of Saunders et al. (1990), 
assuming a dust temperature of 40 K, a dust emissivity index of /3 = 1.2, and pure luminosity evolution of the form L{z) = L(0)(1 + z)'^ 
out to z = 2 (beyond which the luminosity function is simply frozen), for cosmologies (Om = 1.0, Qa = 0.0) and {Hm = 0.3, Ha = 0.7) 
respectively. The lower solid and dashed curves show the number counts predicted by the same models if no luminosity evolution is 
included. The dotted line shows the counts model assumed when creating simulated images. The dot-dash and dot-dot dash assume a 
more complicated, but arguably more realistic luminosity evolution of the form L(z) = L(0)(1 -I- z)^'-^sech^[b ln(l + z) — cjcosh^c , for 
cosmologies {Hm = 1.0, Ha = 0.0) and {Qm = 0.3,f2A = 0.7) respectively, and assuming a dust temperature of 37K and dust emissivity 
index of /3 = 1.2 (Jameson 2000, Small et al. 2002). 
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Figure 12. 2-point angular correlation function for the ELAIS N2 survey field. w{0) 
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fake > 3.50(7 sources respectively, and DD and DR are the number of distinct data pairs and cross data-random pairs within 29 arcsec 
width bins as a function of angular separation 6. The power-law line indicates the correlation function found by Daddi et al. (2000) for 
EROs with R — K > 5 and K < 18.5 as discussed in Sections 9.3 and 9.4. 
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Figure 13. 2-point angular correlation function for the Lockman Hole East survey field. w{8) — _i /r>R\ 
real and fake > 3.50o" sources respectively, and DD and DR are the number of distinct data pairs and cross data-random pairs within 
29 arcsec width bins as a function of angular separation 6. The power-law line indicates the correlation function found by Daddi et al. 
(2000) for EROs with R — K > 5 and K < 18.5 as discussed in Sections 9.3 and 9.4. Note that the vertical dashed line demonstrates the 
cut-off at which there are no more source pairings in the real data. 
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grains radiating at temperature Td- The main uncertainty 
in determining the dust mass Md is the uncertainty in the 
rest-frequency mass absorption coefficient k]^'^^ ■ We have 
adopted the same approach as Hughes, Dunlop & Rawlings 
(f997), using their average value of fcd(800 ^m)= 0.15 ± 
0.09 m'^kg"^, and interpolating to shorter wavelengths by 
assuming kd oc A^^'^. A different choice of fcd(800 ^m) 
would be expected to change the dust mass estimates by a 
factor of up to 2. 

The results are outlined in Tables 9 and 10, the un- 
bracketted quantities assuming an Einstein-de Sitter cosmol- 
ogy, and the values enclosed in brackets adopting Q.m = 0.3, 
JIa = 0.7. SFRs range from several hundred to several thou- 
sand solar masses per year (exceeding even that of the ex- 
treme local starburst Arp220) , the star- forming activity be- 
ing heavily obscured by 10* — 10^ M0 of dust. In fact, for 
Q.M = 0.3, f^A = 0.7 and Ho = 67 kms"^Mpc"^ it can 
be seen that an 850 fim fiux density of 8 mjy corresponds 
closely to a star-formation rate of 1000 Meyr^^. 

The 19 sources brighter than 8 mJy, with S/N > 3.5 
and located in the uniform noise regions of the survey maps, 
account for ~ 10% of the sub-millimetre background ob- 
served by COBE-FIRAS (Puget et al. 1996, Fixsen et al. 
1998, Hauser et al. 1998). Table 11 shows the inferred star 
formation rate density for a variety of redshift bands, and 
for both of our adopted cosmologies. These bright SCUBA 
sources alone imply a high-redshift star formation rate den- 
sity (SFRD) in the range ~ 0.01-0.07 Moyr'^Mpc"^ , com- 
parable to that observed in the optical/UV (Madau et al. 
1996, Steidel et al. 1999). If we assume that the whole of the 
sub-millimtere background may be attributed to starlight 
reprocessed by dust, then this would imply a high-redshift 
SFRD in the range ~ 0.1 - 0.7 Meyr'^Mpc"^, varying only 
by a factor of 2 — 3 on the adopted redshift band for a given 
choice of cosmology. Our results agree very well with Barger, 
Cowie & Richards (2000) who considered the contribution 
to the SFRD of sources brighter than 6 mJy, and determined 
completeness corrected SFRDs ~ 0.1 — 0.4 Moyr^^Mpc^^ in 
the redshift range 1 < 2 < 3, and ~ 0.1-0.7 Moyr'^Mpc"^ 
in the redshift range 3 < z < 6. In contrast to the marked 
decline in the SFRD at 2 > 1 originally implied by op- 
tical/UV observations (Madau et al. 1996), sub-millimetre 
surveys suggest that the SFRD is either steady or gently 
increasing to perhaps as far back as 2 = 5. Current redshift 
constraints (Dunlop 2001a and references within) suggest 
that only 10 — 15% of luminous (> 4 mJy) sub-millimetre 
sources lie at 2 < 2, and that the median redshift of this 
population is Zn^^d — 3. A peak in the SFRD around this 
epoch would not be unexpected given the strong correlation 
between black-hole and spheroid mass found at low redshift 
(Kormendy & Gebhardt 2001) and the peak in optical emis- 
sion from powerful quasars at 2 ~ 2.5. However, improved 
redshift constraints are required to establish when (and in- 
deed if) the comoving SFRD reached a maximum. 

9.3 Co-moving number density of 8 mJy sources 

As explained in the previous Section, although we currently 
lack redshifts for the sub-mm sources uncovered in this sur- 
vey, existing constraints mean that we can be very confi- 
dent that virtually all of them lie at 2 > 1, and indeed 
can be reasonably confident that the majority lie at 2 > 2 



(Dunlop 2001a, Fox et al. 2002). It is thus stiU possible to 
make a meaningful estimate of the co-moving number den- 
sity of bright {Ss5o > 8 mJy) sub-mm sources at high red- 
shift implied by the number counts derived from this sur- 
vey. The results of this calculation, for a variety of assumed 
redshift bands, are summarized in Table 12 (again for our 
two standard assumed alternative cosmologies). From this 
table it can be seen that the co-moving number density of 
dust-enshrouded starburst galaxies with star formation rates 
> 1000 Moyr"^ lies in the range ~ 1 - 10 x 10"^ Mpc"^ 
in very good agreement with Lilly et al. (1999), Barger, 
Cowie & Sanders (1999), and Barger, Cowie & Richards 
(2000). In a similar manner to the comoving SFRD, this 
result depends only weakly (a factor of 2 — 3) on the pre- 
cise choice of assumed redshift band for the sources, for 
a given choice of cosmology. Adopting the now strongly- 
favoured fiat, A-dominated cosmology leads us to conclude 
that the co- moving number density of dust-enshrouded star- 
burst galaxies with star formation rates > 1000 Moyr"^ is 
~ 1 X lO"'"' Mpc"^ 

This is an interesting number. It is over an order of mag- 
nitude smaller than the number density of galaxies brighter 
than L* in the present-day universe, as inferred from the K- 
band luminosity function (Glazebrook et al. 1995; Cowie et 
al. 1995; Gardner et al. 1997; Szokoly et al. 1998; Kochanek 
et al. 2001) , but is an order of magnitude greater than the co- 
moving number density of bright optical QSOs (My < —24) 
at 2 ~ 2 - 3 (Warren, Hewett & Osmer 1995). 

With reference to the present-day .ft'-band luminosity 
function, a co-moving number density of 1 x 10~^ Mpc""^ cor- 
responds to a galaxy luminosity 1 — 1.5 magnitudes brighter 
than L* , or equivalently to galaxies 3 — 4 times more mas- 
sive than an evolved L* galaxy. In this regime the present- 
day galaxy population is completely dominated by massive 
ellipticals (Kochanek et al. 2001). Interestingly the hosts 
of present-day FRII radio sources, and the more luminous 
radio-quiet quasars are also confined to this same high-mass 
regime (Dunlop 2001b, Dunlop et al. 2002). 

Thus, if one wants to attempt to link the high-redshift 
population of very luminous (i.e. SFR ~ 1000 Meyr^^) 
dust-enshrouded starburst galaxies to a low-redshift popula- 
tion, purely on the basis of number-density coincidence, then 
the simplest connection is that the bright SCUBA galaxies 
are the progenitors of present-day massive ellipticals with 
stellar masses ~ 1O^^M0. This is not unreasonable, given 
that such objects require star-formation to be sustained at 
~ 1000 MQyr~^ for ~ 1 Gyr to assemble their present-day 
stellar populations. It is also consistent with the discovery 
that bright SCUBA detections of (comparably massive) ra- 
dio galaxies are largely confined to 2 > 2 (Archibald et al. 
2001). 

If the luminous sub-mm SCUBA sources uncovered by 
this survey are indeed the progenitors of the present-day 
massive elliptical population (Dunlop 2001c), and if, as 
current evidence suggests, the majority of bright sub-mm 
sources lie at 2 > 2 (Dunlop 2001a), then one might rea- 
sonably expect to find a comparably-numerous population 
of passively-evolving massive ellipticals at intermediate red- 
shifts. One way to search for such a population is via surveys 
designed to detect extremely red objects (EROs), such as 
that recently undertaken by Daddi et al. (2000, 2002). One 
can obtain a rough estimate of the co-moving number den- 
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Source 


RA 


DEC 


S'sso 


S/N 


Noise Region 


logioLeo 


SFR 


log 10 


VIdust 




(J2000) 


(J2000) 


/mjy 


850 /im 


850 /im 


/Lq 


/Moyr-i 


/Mq 


01 


16:37:04.3 


41:05:30 


11.2 ±1.6 


8.59 


uniform 


12.49 


12.83) 


649 


1430) 


8.62 


8.96) 


02 


16:36:58.7 


41:05:24 


10.7 ±2.0 


6.27 


uniform 


12.47 


12.82) 


623 


1372) 


8.60 


8.95) 


03 


16:36:58.2 


41:04:42 


8.5 ±1.6 


5.86 


uniform 


12.37 


12.71) 


494 


1087) 


8.50 


8.84) 


04 


16:36:50.1 


40:57:33 


8.2 ±1.7 


5.18 


uniform 


12.36 


12.70) 


480 


1056) 


8.49 


8.83) 


05 


16:36:35.6 


40:55:58 


8.5 ±2.2 


4.16 


uniform 


12.37 


12.71) 


493 


1085) 


8.50 


8.84) 


06 


16:37:04.2 


40:55:45 


9.2 ±2.4 


4.13 


uniform 


12.41 


12.75) 


537 


1182) 


8.54 


8.88) 


07 


16:36:39.4 


40:56:38 


9.0 ±2.4 


4.07 


uniform 


12.40 


12.74) 


523 


1152) 


8.53 


8.87) 


08 


16:36:58.8 


40:57:33 


5.1 ±1.4 


3.82 


uniform 


12.15 


12.49) 


298 


656) 


8.28 


8.62) 


09 


16:36:22.4 


40:57:05 


9.0 ±2.5 


3.76 


non-uniform 


12.40 


12.74) 


526 


1159) 


8.53 


8.87) 


10 


16:36:48.8 


40:55:54 


5.4 ±1.5 


3.69 


uniform 


12.18 


12.52) 


315 


694) 


8.31 


8.65) 


11 


16:36:44.5 


40:58:38 


7.1 ±2.0 


3.67 


uniform 


12.29 


12.63) 


411 


904) 


8.42 


8.76) 


12 


16:37:02.5 


41:01:23 


5.5 ±1.6 


3.65 


uniform 


12.18 


12.52) 


319 


702) 


8.31 


8.65) 


13 


16:36:31.2 


40:55:47 


6.3 ±1.9 


3.56 


uniform 


12.24 


12.58) 


365 


804) 


8.37 


8.71) 


14 


16:36:19.7 


40:56:23 


11.2 ±3.3 


3.55 


non- uniform 


12.49 


12.83) 


651 


1434) 


8.62 


8.96) 


15 


16:37:10.2 


41:00:17 


5.0 ±1.5 


3.52 


uniform 


12.14 


12.48) 


290 


639) 


8.27 


8.61) 


16 


16:36:52.3 


41:05:52 


12.9 ±3.9 


3.51 


non-uniform 


12.55 


12.90) 


751 


1655) 


8.68 


9.03) 


17 


16:36:51.4 


41:05:06 


5.7±1.7 


3.50 


uniform 


12.20 


12.54) 


331 


728) 


8.33 


8.67) 


18 


16:36:11.4 


40:59:26 


20.8 ±6.2 


3.49 


non-uniform 


12.76 


13.10) 


1212 


2669) 


8.89 


9.23) 


19 


16:36:35.9 


41:01:38 


9.4 ±2.8 


3.49 


uniform 


12.42 


12.76) 


550 


1211) 


8.55 


8.89) 


20 


16:36:49.4 


41:04:17 


6.6 ±2.0 


3.48 


uniform 


12.26 


12.60) 


383 


844) 


8.39 


8.73) 


21 


16:37:10.5 


41:00:49 


4.9 ±1.5 


3.47 


uniform 


12.13 


12.47) 


283 


624) 


8.26 


8.60) 


22 


16:36:27.9 


40:54:04 


13.4 ±4.1 


3.46 


non-uniform 


12.57 


12.91) 


780 


1718) 


8.70 


9.04) 


23 


16:37:19.5 


41:01:38 


12.4 ±3.8 


3.45 


non-uniform 


12.54 


12.88) 


721 


1587) 


8.67 


9.01) 


24 


16:36:26.9 


41:02:23 


10.4 ±3.2 


3.43 


non-uniform 


12.46 


12.80) 


607 


1336) 


8.59 


8.93) 


25 


16:36:18.3 


40:59:12 


12.1 ±3.8 


3.37 


non- uniform 


12.52 


12.87) 


703 


1548) 


8.65 


9.00) 


26 


16:36:32.0 


41:00:05 


9.7±3.1 


3.26 


uniform 


12.43 


12.77) 


567 


1248) 


8.56 


8.90) 


27 


16:36:48.3 


41:03:52 


6.6 ±2.1 


3.25 


uniform 


12.26 


12.61) 


386 


851) 


8.39 


8.74) 


28 


16:36:47.2 


41:04:48 


6.3 ±2.0 


3.24 


uniform 


12.24 


12.59) 


367 


809) 


8.37 


8.72) 


29 


16:36:24.1 


40:59:35 


9.9 ±3.3 


3.14 


uniform 


12.44 


12.78) 


573 


1263) 


8.57 


8.91) 


30 


16:37:07.5 


41:02:37 


4.6 ±1.5 


3.13 


uniform 


12.11 


12.45) 


268 


590) 


8.24 


8.58) 


31 


16:36:28.1 


41:01:41 


6.8 ±2.3 


3.07 


uniform 


12.27 


12.62) 


393 


866) 


8.40 


8.75) 


32 


16:36:39.8 


41:00:34 


5.4 ±1.8 


3.07 


uniform 


12.18 


12.52) 


315 


694) 


8.31 


8.65) 


33 


16:36:50.5 


40:58:54 


4.8 ±1.6 


3.06 


uniform 


12.12 


12.47) 


280 


616) 


8.25 


8.60) 


34 


16:36:27.0 


40:58:15 


6.9 ±2.3 


3.05 


uniform 


12.28 


12.62) 


400 


881) 


8.41 


8.75) 


35 


16:36:44.8 


40:56:51 


5.5 ±1.9 


3.02 


uniform 


12.18 


12.52) 


318 


700) 


8.31 


8.65) 


36 


16:36:52.9 


41:02:52 


3.8 ±1.3 


3.00 


uniform 


12.02 


12.37) 


222 


489) 


8.15 


8.50) 



Table 9. Inferred FIR/60 lira luminosities, star-formation rates and dust masses for the sources in ELAIS N2. The un-bracketted 
numbers in columns 7, 8 and 9 assume an Einstein - de Sitter cosmology, whereas the numbers enclosed in brackets assume VIm = 0.3 
and f^A = 0.7. A dust temperature of 40K, and ernissivity index (3 = 1.2 were assumed throughout, as was Hq = 67 krns^^Mpc^-'^. 



sity of passively-evolving massive ellipticals at intermediate 
redshifts by considering the surface density of EROs in the 
Daddi et al. survey with R—K > 5.3 (setting a lower redshift 
boundary of z ~ 1) and K < 18.5 (setting an approximate 
upper redshift boundary of 2: ~ 2 for ellipticals of compa- 
rable mass to bright radio galaxies - Jarvis et al. (2001)). 
The surface density of such objects is ~ 0.1 per sq. arcmin 
(= 350 per sq. degree), very similar to the surface density 
found here for 8 mJy SCUBA sources. For the appropriate 
redshift band 1 < 2: < 2 this surface density converts into 
a co-moving number density of 3 x 10^^ Mpc^"' assuming 
Qm = 0.3, Qa = 0.7 which, as can be seen from Table 12, 
agrees with the co-moving number density of bright SCUBA 
sources to within a factor of 2 or 3 (depending on the choice 
of assumed redshift band for the SCUBA sources) . 

This numerical coincidence provides further circum- 
stantial evidence for an evolutionary path for all mas- 
sive ellipticals which mirrors that which has been already 



largely established for massive radio galaxies (Archibald et 
al. 2000a; Willott et al. 2001; Jimenez et al. 1999), i.e. 



SCUBA source at z > 2.5 ^ 

ERO at 2 ~ 1.5 -* 

3 — 4L* evolved elliptical at z 



A key test of this picture will be to establish whether 
or not the bright SCUBA sources display comparable or 
even stronger spatial clustering than the EROs, for which 
Daddi et al. (2000) report ro — 11 h~^Mpc (for objects with 
R ~ K > 5). As described in the next Section, there is tan- 
talyzing evidence of angular clustering in the 8-mJy survey, 
but it is clear that a substantially larger sub-mm survey 
(approaching 1 degree in size) will be required to settle this 
issue. 
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Source 



RA 
(J2000) 



DEC 
(J2000) 



•Ssso 
/mjy 



S/N 
850 urn 



Noise Region 
850 /im 



logioLeo 

/L0 



SFR 

/Moyr-i 



logloMdus 

/M0 



01 


10:52:01.4 


57:24:43 


10.5 ±1.6 


8.10 


deep strip 


12.47 


12.81) 


613 


1350) 


8.60 


8.94) 


02 


10:52:38.2 


57:24:36 


10.9 ±2.4 


5.22 


uniform 


12.48 


12.82) 


635 


1399) 


8.61 


8.95) 


03 


10:51:58.3 


57:18:01 


7.7 ±1.7 


5.06 


deep strip 


12.33 


12.67) 


448 


987) 


8.46 


8.80) 


04 


10:52:04.1 


57:25:28 


8.3 ±1.8 


5.03 


deep strip 


12.36 


12.71) 


484 


1067) 


8.49 


8.84) 


05 


10:51:59.3 


57:17:18 


8.6 ±2.0 


4.57 


deep strip 


12.38 


12.72) 


500 


1100) 


8.51 


8.85) 


06 


10:52:30.6 


57:22:12 


11.0 ±2.6 


4.50 


uniform 


12.48 


12.83) 


638 


1405) 


8.61 


8.96) 


07 


10:51:51.5 


57:26:35 


8.1 ±1.9 


4.50 


deep strip 


12.35 


12.70) 


474 


1044) 


8.48 


8.83) 


08 


10:52:00.0 


57:24:21 


5.1 ±1.3 


4.38 


deep strip 


12.15 


12.49) 


295 


649) 


8.28 


8.62) 


09 


10:52:22.7 


57:19:32 


12.6 ±3.2 


4.20 


uniform 


12.54 


12.88) 


736 


1620) 


8.67 


9.02) 


10 


10:51:42.4 


57:24:45 


12.2 ±3.1 


4.18 


uniform 


12.53 


12.87) 


712 


1567) 


8.66 


9.00) 


11 


10:51:30.6 


57:20:38 


13.5 ±3.5 


4.09 


non-uniform 


12.57 


12.92) 


786 


1731) 


8.70 


9.05) 


12 


10:52:07.7 


57:19:07 


6.2 ±1.6 


4.01 


deep strip 


12.23 


12.58) 


358 


789) 


8.36 


8.71) 


13 


10:51:33.6 


57:26:41 


9.8 ±2.8 


3.69 


uniform 


12.43 


12.78) 


569 


1252) 


8.56 


8.91) 


14 


10:52:04.3 


57:26:59 


9.5 ±2.8 


3.61 


uniform 


12.42 


12.76) 


552 


1215) 


8.55 


8.89) 


15 


10:52:24.6 


57:21:19 


11.7±3.4 


3.60 


uniform 


12.51 


12.85) 


680 


1497) 


8.64 


8.98) 


16 


10:52:27.1 


57:25:16 


6.1 ±1.8 


3.56 


uniform 


12.23 


12.57) 


356 


784) 


8.36 


8.70) 


17 


10:52:16.8 


57:19:23 


9.2 ±2.7 


3.55 


uniform 


12.41 


12.75) 


537 


1182) 


8.54 


8.88) 


18 


10:51:55.7 


57:23:12 


4.5 ±1.3 


3.55 


deep strip 


12.09 


12.44) 


261 


575) 


8.22 


8.57) 


19 


10:52:29.7 


57:26:19 


5.5 ±1.6 


3.54 


uniform 


12.18 


12.53) 


321 


708) 


8.32 


8.66) 


20 


10:52:37.7 


57:20:30 


10.3 ±3.1 


3.51 


non-uniform 


12.46 


12.80) 


599 


1319) 


8.59 


8.93) 


21 


10:52:01.7 


57:19:16 


4.5 ±1.3 


3.50 


deep strip 


12.10 


12.44) 


262 


577) 


8.22 


8.57) 


22 


10:52:05.7 


57:20:53 


4.7 ±1.4 


3.49 


deep strip 


12.11 


12.46) 


272 


599) 


8.24 


8.59) 


23 


10:51:47.0 


57:24:51 


7.4 ±2.2 


3.48 


uniform 


12.31 


12.66) 


431 


950) 


8.44 


8.79) 


24 


10:51:42.9 


57:24:12 


11.4 ±3.4 


3.47 


uniform 


12.50 


12.84) 


661 


1456) 


8.63 


8.97) 


25 


10:52:36.0 


57:18:20 


12.7 ±3.8 


3.46 


non-uniform 


12.55 


12.89) 


740 


1630) 


8.68 


9.02) 


26 


10:52:27.3 


57:19:06 


8.3 ±2.6 


3.39 


uniform 


12.36 


12.71) 


486 


1070) 


8.49 


8.84) 


27 


10:51:53.8 


57:18:47 


5.5 ±1.7 


3.38 


deep strip 


12.18 


12.52) 


318 


701) 


8.31 


8.65) 


28 


10:52:34.6 


57:20:02 


10.2 ±3.2 


3.31 


uniform 


12.45 


12.79) 


594 


1308) 


8.58 


8.92) 


29 


10:52:16.4 


57:25:07 


6.7±2.1 


3.30 


uniform 


12.27 


12.61) 


389 


856) 


8.40 


8.74) 


30 


10:52:42.2 


57:18:28 


11.2 ±3.6 


3.25 


non-uniform 


12.49 


12.84) 


653 


1437) 


8.62 


8.97) 


31 


10:52:03.9 


57:20:07 


4.0 ±1.3 


3.24 


deep strip 


12.05 


12.39) 


234 


516) 


8.18 


8.52) 


32 


10:52:00.0 


57:20:39 


4.3 ±1.4 


3.22 


deep strip 


12.08 


12.42) 


253 


556) 


8.21 


8.55) 


33 


10:51:33.8 


57:19:29 


8.1 ±2.6 


3.20 


uniform 


12.35 


12.69) 


472 


1040) 


8.48 


8.82) 


34 


10:52:09.9 


57:20:40 


8.5 ±2.8 


3.16 


uniform 


12.37 


12.72) 


495 


1091) 


8.50 


8.85) 


35 


10:51:57.6 


57:26:03 


6.7±2.3 


3.02 


uniform 


12.27 


12.61) 


392 


864) 


8.40 


8.74) 


36 


10:52:03.5 


57:16:54 


8.0 ±2.8 


3.00 


uniform 


12.34 


12.69) 


463 


1019) 


8.47 


8.82) 



Table 10. Inferred FIR/60 ^m luminosities, star-formation rates and dust masses for the sources in the Lockman Hole East. The 
un-brackctted numbers in columns 7, 8 and 9 assume an Einstein - dc Sitter cosmology, whereas the numbers enclosed in brackets 
assume Q,m = 0.3 and f^A = 0.7. A dust temperature of 40K, and cmissivity index (3 = 1.2 were assumed throughout, as was Ho = 
67 kms'^Mpc^^ 



Rcdshift range 



Co-moving SFR density 

(Hm = 1.0, Ha = 0.0) 

/M0yr-lMpc-3 



Co-moving SFR density 

{^M = 0.3, Oa = 0.7) 

/M0yr-lMpc-3 



1-3 
2-3 
2-4 
2-5 
3-4 
3-5 



0.028 
0.057 
0.030 
0.022 
0.065 
0.035 



0.018 
0.033 
0.017 
0.012 
0.035 
0.018 



Table 11. The co-moving star formation rate density of sources brighter than 8 mJy, detected with S/N > 3.50, using the uniform noise 
regions of the two survey areas and assuming that all of these sources lie within the redshift range given in column 1 . Column 2 assumes 
an Einstein-de Sitter cosmology, and column 3 assumes Qm = 0.3, S^a = 0.7. Hq = 67 kms~^Mpc~^ was adopted for both cosmologies. 
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Rcdshift range 



Co-moving no. density 

(Q.M = 1-0, Ha = 0.0) 
/Mpc-3 



Co-moving no. density 

(Qm = 0.3, Qa = 0.7) 

/Mpc-3 



1-3 


4.94 X IQ- 


1.40 X 10--' 


2-3 


9.96 X IQ- 


2.63 X 10--'' 


2-4 


6.05 X 10- 


1.34 X 10-"' 


2-5 


3.74 X lo- 


0.94 X 10-"' 


3-4 


ll. 22 X 10- 


2.77 X 10-"' 


3-5 


6.00 X 10- 


1.45 X 10-"' 



Table 12. The co-moving number density of sources brighter than 8 mjy, detected with S/N > 3.50, using the uniform noise regions 
of the two survey areas and assuming that all of these sources lie within the redshift range given in column 1. Column 2 assumes an 
Einstein-de Sitter cosmology, and column 3 assumes Qm = 0.3, Qa = 0.7. Hq = 67 kms-^Mpc-^ was adopted for both cosmologies. 



9.4 Clustering 

If the bright 850 /im sources are indeed the progenitors 
of massive elliptical galaxies then they should be strongly 
clustered, an inevitable result of gravitational collapse from 
Gaussian initial density fluctuations since the rare high-mass 
peaks are strongly biased with respect to the mass. There is 
a great deal of evidence to support the presence of this bias 
at high redshift. The correlations of Lyman-break galaxies 
at z ~ 3 (Steidel et al. 1999) are almost identical to to those 
of present-day field galaxies, even though the mass must 
have been much more uniform at early times. Furthermore, 
the correlations increase with UV luminosity (Giavalisco & 
Dickinson 2001) reaching scale lengths of r^ ~ 7.5h~^ Mpc 
- approximately 1.5 times the present-day value. In the case 
of luminous proto-ellipticals we expect an even stronger bias 
since we select not just massive galaxies but those that have 
collapsed particularly early in order to generate the oldest 
stellar populations. This is suggested by studies of the lo- 
cal Universe which have shown that early-type galaxies are 
much more clustered than late- type galaxies (eg. Guzzo et al. 
1997, Willmer et al. 1998), and more recently by the findings 
of Daddi et al. (2000) who have investigated the clustering 
properties of extremely red objects (EROs). They detect a 
strong clustering signal of the EROs which is about an or- 
der of magnitude larger than the clustering of if— selected 
field galaxies, and also report a smooth trend of increasing 
clustering amplitude with increasing R — K colour, reaching 
ro — ll/i~^ Mpc for R — K > 5. These results are probably 
the strongest evidence to date that the largest fraction of 
EROs is composed of ellipticals at 2 > 1. 

We already have some hints of the possibility of strong 
clustering in the bright sub-mm population from the dis- 
covery of a strong excess of bright SCUBA sources around 
high-redshift AGN (Ivison et al. 2000). Motivated by this 
result, and the potential connection with the intermediate- 
redshift EROs discussed above, we have attempted to quan- 
tify the possible strength of clustering in our two survey 
fields via calculation of the the two-point correlation func- 
tion. The results, based purely on sources with S/N > 3.50, 
are shown in Figs 12 and 13. There is tentative evidence 
of clustering on scales of 1 — 2 arcmin in both of our sur- 
vey fields, but most particularly, and somewhat stronger, in 
the ELAIS N2 region. Referring to Fig 9, we can in fact see 
by eye that the most significant {S/N > 3.50) sources in 
ELAIS N2 do not conform to a homogeneous distribution 
across the field - rather there are two apparent concentra- 



tions of 850 /im sources in the top left and bottom right of 
the image, approximately at RA 16:37:00, DEC -f41:05:00, 
and RA 16:36:30, DEC -f40:56:00 respectively, with an ap- 
parent under-density of sources in the intervening regions. A 
second peak in ■w{6) at ~ 600 arcsec (the distance between 
these over-densities), and a trough at ~ 400 arcsec are seen 
to reflect this "by eye" distribution. Very interestingly, we 
see the same large-scale inhomogeneities in the Chandra X- 
ray image of the ELAIS N2 region (Almaini et al. 2002), 
although the coincidence of X-ray and SCUBA sources is 
small (< 10%). Almaini et al. provide a full discussion and 
propose a possible explanation for this effect. 

It is clear from Figs 12 and 13 that, due to small-number 
statistics, this first attempt at a direct measure of the clus- 
tering of bright sub-mm sources has proved inconclusive. 
However, this should not be taken as evidence that the 
SCUBA sources are unclustered. In fact, as illustrated in 
these figures, it is worth noting that these correlation func- 
tions are certainly still consistent with the strong clustering 
signal detected for EROs by Daddi et al. (2000). Clearly a 
much larger survey (~ 0.5 square degrees) will be required to 
obtain a meaningful measurement of the strength of cluster- 
ing in the bright sub-mm population. Near complete redshift 
information may also be required to quantify the extent to 
which any clustering signal will be partially erased by projec- 
tion through a wide range in redshift. For example, suppose 
that the bright SCUBA population spans a relatively wide 
range in redshift 2 < 2 < 5 whereas, as argued above, the 
bright EROs are confined to the redshift range 1 < 2 < 2. In 
that case, if the strength of the spatial clustering in the two 
populations was the same, w{6) as measured from SCUBA 
images would be expected to be 1.6 — 1.9 times smaller (de- 
pending on choice of cosmology) than that which has been 
measured for EROs. 



10 CONCLUSIONS 

We have used the SCUBA camera on the JCMT to con- 
duct the largest extra-galactic sub-mm survey to date, in 
two regions of sky (ELAIS N2 and Lockman Hole East) 
covering a total area of 260 arcmin^ to a typical rms noise 
level of (7850 — 2.5 mjy/beam. We have reduced the re- 
sulting dataset using both the standard JCMT SURF soft- 
ware and our own IDL-based pipeline which allows us to 
produce zero-footprint maps and noise images. The uncor- 
related noise maps from the latter approach have enabled 
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us to quantify the statistical significance of each peak in our 
maps leading to properly quantified errors on the flux den- 
sity of our sources. We find f 9 sources with S/N > 4.00, 38 
sources with S/N > 3.50 and 72 sources with S/N > 3.00. 
A series of Monte-Carlo simulations have allowed us to as- 
sess the completeness, positional uncertainties and effects 
of confusion in this survey, leading to the most accurate 
estimate of the sub-millimetre counts over the flux-density 
range 5*850 = 5 ~ f2 mjy. Our best estimate of the cumu- 
lative source count at Sgso > 8 mJy is 320^iqq, consistent 
with a range of models involving strong evolution of the 
dust-enshrouded starburst population. 

Adopting Q,M = 0.3, Qa = 0.7 and Hq — 
67 kms~^Mpc^^, the sources uncovered by this survey have 
inferred star-formation rates in excess of lOOOMoyr^^, en- 
shrouded in 10* — 10® M0 of dust. Assuming simply that the 
majority of these sources lie at z > 1.5, this result implies 
that their co- moving number density is ~ 10~^Mpc~'', co- 
incident with the number density of massive ellipticals with 
L > 3 — 4L* in the present-day Universe. This co-moving 
number density is also the same as that inferred for compara- 
bly massive passively-evolving objects in the redshift band 
1 < z < 2 from recent surveys of EROs. Thus while the 
bright submillimetre sources uncovered by this survey con- 
tribute only ~ 10% of the submillimetre background, they 
can plausibly account for the formation of all present-day 
massive ellipticals. A key test of this picture is to determine 
whether the bright SCUBA population is at least as strongly 
clustered as are the EROs at intermediate redshifts. As a 
first attempt at such a test we have compared the distribu- 
tion of our sources with a mock catalogue of fake sources. 
Our results have proved inconclusive due to the small num- 
ber of significant sources (S/N> 3.50) detected in our survey 
areas, but there are some initial indications of clustering on 
scales of 1 — 2 arcminutes, particularly in the ELAfS N2 
field. 

The process of multi- wavelength follow-up of this survey 
is already well underway. SED-based redshift constraints, 
along with potential optical-infrared identifications of our 
most significant sources are presented in paper II (Fox et al. 
2002), while the results of IRAM PdB interferometer imag- 
ing of the most significant Lockman Hole East is presented 
by Lutz et al. (2001). Correlations with deep X-ray and ra- 
dio images will also be discussed in subsequent papers by 
Almaini et al. (2002) and Ivison et al. (2002). 
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APPENDIX A: THE SOURCE-EXTRACTION ALGORITHM. 



Suppose we consider a normalised beam-map B{x,y) as a source template and that at position {i,j) in the unconvolved 
image the signal is S{i,j) and the noise is N{i,j). If n peaks above a specified flux threshold are located in the Gaussian- 
convolved image, we may construct a model to the unconvolved image such that beam-maps centred on each peak position 
are simultaneously scaled to give an overall best fit to the entire image. Using a minimised x^ fit as the maximum likelihood 
estimator then 

2_ [yr^ S{i,j)-Y,l^-^akBk{x~i,y~j)\ 
where at is the best fit fiux to the fcth peak. Minimising with respect to each at 

and m = 1 • • • n. 

Rearranging, we obtain the matrix equation 

E/ V^ Bk[x - i,y - j)Bm(x - i,y - j) \ _ \^ S{i,j)Bm[x-i,y~j) 
which may be written in the form 

n 
y^ amfcOfe = Pm (A4) 

fc=l 
where 



_^ Bk{x-i,y-j)Bm{x~i,y~j) 



E ^k(x-i,y- j)iim(x ~i,y~j} 



an n X n matrix, and 



« _^ S(i,j)B^{x^i,y-j) 



a vector of length n. 

To recover the best fit values of at we invert matrix a such that 



flfc 



^Hfc^/3m (A7) 



The variance associated with the estimate a^ is given by 

Since a^rn is independent of S(i,j) 
dak v^ [a]k^Brr,{x-i,y-j 



r 1 — 1 7-1 / \ 

_ V Mfc„gm(a^-»,t;-j) , , 

dS{i,j) 2^ N{i,jY ^ ^^ 

m — 1 

Consequently we find, 

fe=i i=i \ i,j ^ '■^' ) 

The final term in brackets is simply the matrix [a] and so this expression reduces to 

N\ak) = H,-,^ (All) 

The diagonal elements of [q] ^ are the variances of the fitted parameters a^ such that the significance of the peak detection 
a(a,) = ^= (A12) 

vMfefc 



© 0000 RAS, MNRAS 000, 000-000 



